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ABSTRACT 


Demography  and  genetics  of  increasing  populations  of  blue  grouse 
(Dendragapus  obsourus)  were  studied  by  combining  intensive  field  observations 
with  analysis  of  genetic  structure  of  the  populations.  This  research  was 
designed  to  test  the  hypothesis  that  genetic  structure  of  populations  changes 
with  changes  in  population  density. 

Blue  grouse  colonize  newly  logged  regions  rapidly  and,  in  this  study, 
density  of  breeding  birds  increased  for  at  least  seven  years  following 
logging.  The  age-structure  of  colonizing  populations  changed  from  70%  to 
45%  yearlings  in  the  first  three  years  after  logging  and  stabilized  there¬ 
after. 

The  Ng  locus  appears  as  white  bands  on  gels  stained  for  serum  esterase, 
following  electrophoresis.  The  Ng  locus  was  polymorphic  and  was  used  as  a 
marker.  The  three  alleles  at  this  locus  were  widespread  spatially  and  stable 
temporally.  Frequency  of  the  most  common  allele,  NgM3  was  higher  in  birds 
in  one-year-old  habitat  than  in  birds  on  older-aged  habitats.  Frequency  of 
heterozygotes  as  well  as  the  observed  frequency  of  heterozygotes/expected 
frequency  (0/E  ratio)  among  yearlings  were  significantly  positively  regressed 
on  age  of  habitat.  Correlation  between  population  density  and  both  frequency 
of  heterozygotes  and  the  0/E  ratio  among  yearlings  was  positive  and  signifi¬ 
cant,  i.e.,  the  higher  the  population  density,  the  higher  the  frequency  of 
both  heterozygotes  and  the  0/E  ratio  among  yearlings.  These  results  support 
the  original  hypothesis  that  changes  in  genetic  structure  occur  with  changes 
in  density.  The  cause  of  these  relationships  is  not  known  but  experiments 
were  suggested  to  investigate  two  possible  alternatives.  A  conceptual  model, 
based  on  the  idea  of  total  heterozygosity  of  individuals  conferring  competi¬ 
tive  advantage,  was  described  to  account  for  the  observed  results. 
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"No  population  ecology  can  succeed  that  is  not  also  population 


genetics . " 


Richard  C.  Lewontin  (1968)  in 
Population  Biology  and  Evolution 
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INTRODUCTION 

Populations  of  blue  grouse  on  Vancouver  Island  increase  rapidly 
following  logging  or  burning  of  mature  forest  (Hatter  1955)  suggesting  a 
surplus  of  grouse  is  available  for  colonization.  If  rapid  expansions  are 
possible  in  these  situations,  what  are  the  mechanisms  which  stop  these 
increases?  Chitty  (1965,  1967)  proposed  that  changes  in  the  genetic  quality 
of  populations  associated  with  changes  in  abundance  were  both  necessary  and 
sufficient  to  stop  increases.  Much  recent  research  in  population  biology 
has  attempted  to  see  if  genetic  changes  are  associated  with  changes  in 
abundance. 

The  objectives  of  the  research  reported  in  this  thesis  were  twofold: 

a)  to  describe  the  mechanism  of  population  increase  following  logging,  and 

b)  to  test  the  hypothesis  that  genetic  changes  are  associated  with  changes 
in  abundance.  (The  null  hypothesis  is  that  genetic  changes  are  not 
associated  with  changes  in  density.)  Demography  was  studied  via  intensive 
field  observation  and  genetics  were  studied  using  electrophoresis  to  dif¬ 
ferentiate  genetic  variants  at  a  marker  locus.  An  experimental  approach 
was  made  possible  because  logging  of  mature  forest  continually  creates  new 


habitat. 
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BACKGROUND 


The  Experimental  Approach 

The  creation  of  new,  previously  unoccupied,  breeding  habitat  may  pre¬ 
sent  ecological  opportunities  for  certain  species  and  can  be  used  as  an 
experimental  approach  to  induce  population  increases.  Often,  new  breeding 
habitat  is  an  unexpected  by-product  of  man’s  alteration  of  the  ecosystem 
(Bustard  1969,  Tevis  1956a, b).  Logging  on  Vancouver  Island,  for  example, 
induces  rapid  population  increases  in  blue  grouse  by  creating  favorable 
breeding  habitat.  Present  day  logging  patterns  on  Vancouver  Island  involve 
the  complete  removal,  i.e.  clear- cutting ,  of  blocks  of  forest.  Often 
adjacent  blocks  of  forest  are  cut  in  successive  years  creating  a  mosaic  of 
habitat  types  within  a  restricted  region.  This  experimental  manipulation 
is  not  the  only  way  to  create  new  breeding  habitat  but  it  is  the  most 
practical  and  probably  simulates  fire  in  earlier  times,  at  least  to  grouse. 

The  approach  used  was  to  choose  three  large  areas  containing  sub-units 
which  were  logged  in  several  previous  years.  I  will,  throughout  this  thesis, 
speak  of  age  of  habitat  synonymously  with  years  since  logging.  The  number 
of  grouse  on  these  areas  was  ascertained  and  habitat-specific  densities  were 
calculated  by  associating  each  grouse  with  the  habitat  type.  All  areas 
logged  in  the  same  year  were  combined  for  analysis. 

The  Animal 

Blue  grouse  are  endemic  to  mountains  of  western  North  America.  On 
Vancouver  Island,  in  spring  and  summer,  they  occur  in  openings  in  the  forest 
and  in  winter  they  occur  in  mature  or  old  second-growth  forest  (Bendell  and 
Elliott  1967).  In  spring,  adult  males  (2-year-olds  and  older)  are  territor¬ 
ial  and  live  solitarily.  Yearling  males  (1  year  old)  are  not  usually 
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territorial,  roam  over  extensive  areas  in  spring  and  probably  do  not  contri¬ 
bute  to  breeding.  Females,  both  adults  and  yearlings,  do  not  appear  to  be 
territorial  in  spring  or  summer.  Both  adult  and  yearling  females  breed. 

Most  breeding  occurs  in  May  and  chicks  usually  hatch  in  June.  Grouse  appear 
polygamous  and  males  do  not  participate  in  nesting,  incubation  or  care  of 
the  young.  Once  chicks  hatch,  females  and  broods  move  rather  extensively. 

Grouse  are  suitable  for  study  since  they  are  diurnal,  conspicuous  and 
easily  captured,  banded  and  resighted.  One  drawback  associated  with  studies 
on  blue  grouse  is  their  migration  off  breeding  range  into  dense  forest  in 
late  summer  and  autumn  which  makes  it  difficult  to  study  them  in  winter 
(King  1971).  Other  problems  involve  overlapping  generations,  secretive 
nature  of  1-year-old  males  and  difficulty  in  finding  nests.  However,  these 
are  of  minor  hindrance  to  this  study. 

Research  on  population  dynamics  of  grouse  on  logged  regions  of  Vancouver 
Island  has  been  relatively  intense  since  the  early  1950’s.  Bendell  (1955) 
studied  grouse  near  Campbell  River  and  concluded  that  parasites  were  limiting 
his  very  dense,  but  stable,  population.  In  later  s tudies  he  was  unable  to 
confirm  this.  A  new  study,  started  in  1958,  was  aimed  at  seeing  if  surplus 
grouse  were  available  to  colonize  depopulated  areas.  These  new  studies  were 
on  a  stable  population,  but  densities  were  less  than  half  those  found  in  the 
earlier  study.  Basically,  from  this  study  Bendell  concluded  that  there  were 
no  surplus  adult  grouse  (Bendell  and  Elliott  1967). 

Zwickel  (1965)  and  Zwickel  and  Bendell  (1967a)  studied  the  relationship 
between  nutrition,  chick  survival  and  subsequent  population  density.  They 
concluded  that  survival  of  chicks  was  predetermined  in  the  mother  but  that 
the  level  of  production  in  any  one  year  was  more  than  adequate  to  explain 
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recruitment  the  following  year.  Like  the  earlier  studies,  populations  in 
this  study  were  stable  but  at  a  lower  density. 

Thus,  population  densities  of  grouse  appeared  stable,  indicating  they 
were  regulated,  even  though  no  substantial  surplus  of  potentially  breeding 
grouse  could  be  found.  But  surpluses  are  also  indicators  of  regulated 
populations  (Brown  1969,  Watson  and  Moss  1970).  Reproduction  varied  between 
years  whereas  recruitment  did  not,  suggesting  recruitment  was  regulated, 
at  least  in  stable  populations.  Food  seemed  adequate  and  disease  was  un¬ 
common.  Most  external  limiting  factors  seemed  inadequate,  either  by  them¬ 
selves  or  in  combination,  to  explain  numbers  of  grouse  (Zwickel  and  Bendell 
1967a,  1970).  That  leaves,  as  one  possibility,  the  interactions  of  the 
grouse  themselves  as  a  regulatory  mechanism. 

Chitty  (1967)  has  pointed  out  the  problem  of  discovering  a  regulatory 
mechanism  in  a  stationary  population.  Since  stationary  populations  can  offer 
no  comparative  information  on  the  causes  of  increases  and  decreases,  they 
are  of  limited  value  in  the  study  of  population  regulation.  That  is  why 
the  experimental  approach  used  in  this  study  seemed  necessary  in  an  attempt 
to  explain  population  regulation  in  blue  grouse. 

The  Genetic  Marker 

Studies  involving  the  interaction  of  genetics  with  ecology  require  a 
suitable  genetic  marker (s).  Criteria  most  often  used  in  choosing  markers 
are  that  they  be  a)  variable,  b)  easily  identified,  c)  simply  inherited,  and 
d)  stable.  In  the  past,  studies  of  ecological  genetics  such  as  this  one  were 
not  possible  with  most  species  since  these  criteria  could  not  be  met.  Even 
though  morphological  variants  were  known  to  occur  in  some  species  [e.g., 
banding  in  Cepaea  memoralis  (Cain  and  Sheppard  1954) ;  melanism  in  Biston 
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betulaP'ia  (Kettlewell  1956);  spot  patterns  in  Maniola  jurtina  (Ford  1965a); 
jP-1ocus  in  Mus  musculus  (Lewontin  and  Dunn  I960)],  it  was  not  until  the 
development  of  electrophoresis  that  a  tremendous  amount  of  cryptic  genetic 
variation  was  revealed  in  virtually  all  species  (Lush  1966,  1970,  Manwell 
and  Baker  1970) .  When  I  initiated  this  research  nothing  was  known  about 
genetic  variation  within  populations  of  blue  grouse.  However,  several 
proteins  and  enzymes  of  other  vertebrates  were  characteristically  polymorphic 
(Figure  1;  Lush  1967,  1970,  Manwell  and  Baker  1970). 

My  aim,  at  first,  was  to  survey  several  proteins  and  enzymes  to  see 
which  showed  electrophoretic  variation.  Of  those  surveyed,  transferins  and 
lactate  dehydrogenase  were  not  polymorphic  while  some  others  (including  acid 
phosphatase,  alkaline  phosphatase  and  leucine  amino  peptidase)  gave  consis¬ 
tently  poor  resolution.  Esterases  were  variable  but  the  genetics  of  this 
variation  seemed  too  complex.  Baker  et  al,  (1966)  reported  a  similar  com¬ 
plexity  in  the  esterases  of  pheasants  (Phasianus  eolchicus)  in  North  America. 
Since  I  was  searching  for  genetic  variation  which  could  be  investigated 
without  an  elaborate  breeding  program,  I  disregarded  esterases. 

However,  during  a  routine  stain  for  serum  esterase,  early  in  this  study, 
a  gel  was  accidentally  overstained,  and  when  washed  in  the  normal  manner 
(see  Methods) ,  revealed  a  series  of  white  bands  (Figure  2)  on  a  blue  or 
grey  background.  These  bands  ( Ng  locus)  were  variable  and  of  apparent 
simple  inheritance  (Birdsall  et  al ,  1970)  and  were  chosen  for  study  on  this 
basis.  Thus,  the  marker  was  chosen  more  or  less  at  random.  The  biochemical 
significance  of  the  substance  is  unknown.  However,  this  matters  little  to 
the  analysis,  since  the  hypothesis  was  stated  in  general  terms  to  cover  any 
genetic  variation. 
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Figure  1. 


Number  of  known  biochemical  loci  in  vertebrates.  Line  A 
represents  those  loci  discovered  by  the  techniques  of 
electrophoresis;  line  B  represents  loci  discovered  by  other 
techniques  [modified  after  Lush  (1370)]. 
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Starch-gel  showing  the  three  white  bands  and  four 
genotypes  at  the  Ng  locus  in  blue  grouse. 

1.  NgM/NgM ;  2.  NgS/NgM;  3.  NgS/NgM ;  4.  NgS/NgM; 
5.  NgS/NgF ;  6.  NgM/NgF;  7.  NgS/NgM;  8.  NgM/NgM. 


Figure  2. 
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Genetics  and  Demography 

There  has  been  an  increasing  interest  for  the  past  several  years  in 
considering  demography  and  genetics  simultaneously  (Birch  1960,  Cavalli- 
Sforza  1962,  Cavalli-Sforza  and  Bodmer  1971,  Charlesworth  1971,  Chitty  1960, 
Crow  and  Kimura  1970,  Ford  1965a, b.  King  and  Anderson  1971,  Kojima  1971, 
Lewontin  1970,  MacArthur  and  Wilson  1967,  Roughgarten  1971,  Wallace  1968a, b, 
1970,  and  many  others).  This  interest  has  been  especially  spurred  by  recent 
technical  advances  in  the  use  of  electrophoresis  which  makes  the  detection 
and  differentiation  of  genetic  variation  in  most  species  relatively  easy 
(Hunter  and  Markert  1957,  Poulik  1957,  Smithies  1955). 

Among  numerous  recent  hypotheses  (e.g.  Ford  1965a, b,  Howard  1960,  etc.) 
involving  the  interaction  of  genetics  and  demography,  Chitty' s  (1960,  1964, 
1965,  1967,  1970)  is  one  of  the  more  interesting  and  has  some  evidence  in 
favor  of  it.  For  example.  Ford  (1965a)  showed  that  during  periods  of  rapid 
population  increase  in  butterflies  (Lepidoptera)  selection  pressure  was 
relaxed  and  resulted  in  a  burst  of  genetic  variability.  This  was  countered 
by  increased  selection  during  declines.  Wellington  (1960)  demonstrated 
shifts  in  proportions  of  different  types  of  tent  caterpillars  ( Malacosoma 
pluviale)  during  changes  in  abundance.  Baltensweiler  (1968)  and  Clark  et  at, 
(1967  p.124-136)  argued  for  qualitative  shifts  in  populations  of  the  Larch 
tortix  ( Zeivaphera  grisearia)  with  changes  in  abundance.  Although  Wellington 
and  Baltensweiler  did  not  demonstrate  the  genetics  of  the  types  examined, 
it  seems  unlikely  that  the  changes  in  types  were  totally  induced  by  the 
environment,  since  they  existed  at  all  population  densities  examined. 

Perhaps  the  most  outstanding  works  on  demography  and  genetics  of  verte¬ 
brates  are  those  being  done  on  several  species  of  small  rodents  in  North 
America  (Canham  1969,  Gaines  and  Krebs  1971,  Gaines,  Myers  and  Krebs  1971, 
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Myers  and  Krebs  1971,  Rasmussen  1964,  1970,  Sc lander  1970a, b.  Sc lander  and 
Yang  1969,  Tamar in  and  Krebs  1969).  Rasmussen  and  Selander  aimed  their 
research  at  the  genetic  structure  of  populations  and  how  genetic  structure 
was  affected  by  demography  and  behaviour,  while  the  other  studies  cited 
above  attempted  to  study  demography  and  behaviour  as  affected  by  genetics. 
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STUDY  AREAS 

In  1968  a  brief  survey  was  conducted  on  three  regions  and  more  inten¬ 
sive  work  was  confined  to  two  other  main  regions.  Since  1968,  grouse  on 
the  latter  two  regions,  plus  one  new  area,  were  studied.  Since  the  only 
information  used  from  the  three  general  regions  studied  in  1968  is  allelic 
frequencies,  further  descriptions  of  these  plots  are  not  necessary .  Gener¬ 
ally,  all  areas  were  in  a  relatively  young  stage  of  vegetative  regeneration 
following  logging  or  burning.  The  intensive  study  areas  are  described  below. 
All  are  in  the  Ash  River  Valley  near  Elsie  Lake,  midway  between  Courtenay 
and  Port  Albemi,  Vancouver  Island  (Figure  3).  Sizes  and  logging  histories 
of  each  study  area  are  given  in  Table  1. 

Area  104  is  on  a  steep  southwest-facing  slope  ranging  in  elevation 
from  270  to  1070  meters.  Most  of  it  was  logged  since  1962  but  only  part  of 
the  debris  left  after  logging  was  burned.  This  area  probably  has  the  best 
soil  (J.  Holm,  Chief  Forester,  MacMillan- Bloedel  Sproat  Lake  Division,  pers. 
comm . )  of  all  the  intensively  studied  tracts. 

Area  107  is  the  largest  of  the  study  plots  and  its  size  has  continu¬ 
ously  increased  since  1968,  through  the  addition  of  newly  logged  regions 
adjacent  to  the  study  area.  This  region  is  characterized  by  a  series  of 
parallel  rocky  ridges  running  east  and  west  with  elevations  from  200  to  600 
meters. 

Area  108e  was  added  to  the  study  in  1969  and  is  the  smallest  study 
area.  It  has  the  most  uniform  logging  history  and  the  least  relief,  but 
is  similar  in  topography  to  107. 

The  Ash  River  Valley  would  be  classed  in  the  Coastal  Western  Hemlock 
Zone  (Krajina  1969).  Most  of  the  natural  forest  in  the  Ash  River  Valley  was 
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Figure  3. 


Map  of  the  Ash  River  Valley,  Vancouver  Island,  British 
Columbia,  showing  the  locations  of  the  three  main  study 


areas . 
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Table  1. 


History  of  logging  of 
compiled  from  maps  of 
Bloedel,  Ltd.,  Sproat 


the  study  areas  (in  km2), 
logging  patterns  supplied 
Lake  Division. 


Data  were 
by  MacMillan 


Year 

Area 

104 

107 

108e 

pre- 1961 

0.943 

0.826 

m m 

1961 

- 

0.198 

- 

1962 

0.413 

- 

- 

1963 

- 

- 

- 

1964 

0.287 

0.429 

0.457 

1965 

0.376 

0.696 

- 

1966 

0.405 

0.413 

0.583 

1967 

0.324 

0.304 

0.134 

1968 

0.708 

0.494 

- 

1969 

- 

0.822 

- 

1970 

_ 

0.267 

_ 
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removed  (Figure  4).  The  logging  history  of  the  valley,  as  compiled  from 
maps,  is  presented  in  Table  2.  In  recent  years,  15  to  40  hectare  stands 
of  trees  were  logged  using  clear-cut  methods,  creating  a  patchwork  of 
various  habitat  types.  It  was  not  unusual,  for  example,  to  have  20- ,  10- , 
5-  and  1-year-old  areas  immediately  adjacent  to  each  other.  Since  1960, 
logging  activities  increased  in  the  valley,  especially  when  compared  to  the 
decade  starting  with  1950.  Almost  all  of  the  logged  region  was  replanted 
with  Douglas  fir  and  these  plantations  were  in  various  stages  of  regenera¬ 
tion. 

Although  vegetation  was  not  studied  in  detail  some  general  statements 
can  be  made.  Fireweed  (Epilobium  angusti folium)  was  probably  the  most 
prominent  herbaceous  plant  in  summer.  Douglas  fir  (Pseudotsuga  menziesii) 
and  western  hemlock  (Tsuga  heterophylla)  were  the  predominant  conifers. 
Judging  from  published  results  of  other  studies  on  Vancouver  Island  (Bendell 
and  Elliott  1967,  Zwickel  and  Bendell  1967),  vegetation  on  the  areas  I 
studied  was  more  luxuriant  with  less  willow  (Salix  spp.j. 

Weather  in  the  Ash  River  Valley  is  influenced  by  the  Beaufort  Mountains 
to  the  east.  Rainfall  is  over  250  cm  annually  with  snow  common  in  winter. 
Summers  are  pleasant  with  hottest  days  over  40  C  and  coolest  nights  near  0  C 
(as  measured  in  weather  screens  at  ground  level).  Long  periods  of  very  hot 
and  dry  weather  are  frequent  between  May  1  and  September  1. 
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Figure  4.  Photograph  of  the  Ash  River  Valley  looking  east  towards 
the  Beaufort  Mountains.  Trees  in  most  of  the  valley 
bottom  were  removed  prior  to  1945  while  the  more  open 
areas  on  the  upper  slopes  have  been  logged  since  1960. 
Note  the  clear-cut  logging  pattern.  Crest  of  these 
mountains  is  near  1800  meters. 
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History  of  logging  in  the  Ash  River  Valley.  Data  were  compiled 
from  maps  of  logging  patterns  provided  by  MacMillan  Bloedel, 
Ltd. ,  Sproat  Lake  Division. 
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METHODS 

Finding  and  Marking  Grouse 

Field  work  began  in  late  April  in  1968  and  in  late  March  or  early  April 
in  other  years.  In  all  years  field  work  terminated  near  September  1.  A 
systematic  search  with  trained  dogs  was  conducted  daily.  Grouse  were  located 
by  one  of  three  methods:  1)  the  dogs  pointed  a  majority  of  the  silent  birds, 
2)  in  spring  and  early  summer  males  on  territories  were  conspicuous  and  often 
located  because  of  their  displays  and  calls,  and  3)  once  eggs  hatched,  hens 
with  chicks  were  often  located  by  their  response  to  an  artificial  chick  call 
given  by  the  observer. 

An  attempt  was  made  to  capture  and  band  all  unmarked  grouse.  Grouse 
were  captured  with  a  7-meter  telescoping  noose  pole  (Zwickel  and  Bendell 
1967b).  Birds  in  the  hand  were  weighed,  bled,  marked  with  unique  color 
combinations  of  leg  bands,  and  released.  Chicks  less  than  25  days  of  age 
were  too  small  to  be  marked  with  leg  bands.  Most  of  the  information  in  this 
thesis  comes  from  marked  grouse.  For  purposes  of  analysis,  the  point  of 
first  sighting  in  a  given  year  was  used  to  assign  a  bird  with  a  habitat  of 
known  age. 

Estimating  Density 

Adult  males  were  relatively  easy  to  count  in  spring  and  early  summer 
since  they  occupy  territories  from  which  they  display  and  mate  (Bendell  and 
Elliott  1967).  Most  territorial  males  were  captured  and  marked  but  some 
were  difficult  to  capture  because  they  were  wild  or  usually  located  in  tall 
trees.  If  unmarked  males  with  characteristic  behaviour  were  consistently 
found  in  a  specific  region,  I  included  them  in  the  estimate  of  density. 

Yearling  males,  unlike  adults,  were  more  difficult  to  count  since  they 
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were  usually  not  territorial  and  did  not  display.  Thus,  I  used  an  indirect 
method  of  analysis  based  on  the  sex-ratio  of  yearling  males  and  yearling 
females  counted  prior  to  the  onset  of  hatching  of  eggs  (a  period  when  this 
ratio  is  least  likely  to  be  biased).  (If  this  ratio  was  1:1,  I  assumed  that 
the  numbers  of  yearling  males  and  females  were  equal.)  During  the  four  years 
of  this  study,  142  yearling  females  and  121  yearling  males  were  captured  in 
this  period  (1  •  18-¥ :  Id*) .  This  ratio  was  not  different  from  an  expected  1:1 
ratio  (x2(1)=1.67;  p=0.1)  and  I  conclude  that  there  were  the  same  number  of 
yearling  males  as  yearling  females.  Thus,  density  of  yearling  males  was 
considered  the  same  as  for  yearling  females. 

Like  most  yearling  males,  females  had  few  displays  and  activities  in 
spring  which  made  them  conspicuous.  Likewise,  they  did  not  appear  to  be 
territorial.  However,  once  chicks  hatched,  females  with  broods  were  easy  to 
count  and  most  brood  hens  were  probably  found  in  all  years  (except  perhaps 
1968).  But  not  all  females  had  broods  which  made  females  more  difficult  to 
count.  Zwickel  and  Bendell  (1967a)  used  number  of  females  seen  per  hour  of 
search  during  the  period  May  14  to  June  3  as  an  index  of  density,  but  noted 
that  this  was  potentially  biased.  I  developed  a  method,  modified  from  the 
common  Lincoln  model  (Lincoln  1930) ,  for  estimating  density  which  circumvented 
most  of  the  biases  of  the  method  of  Zwickel  and  Bendell  (1967a). 

Estimation  formulae.  Numbers  of  adult  and  yearling  females  were  esti¬ 
mated  separately,  Na  and  Ny.:  For  purposes  of  illustration,  consider  only 
adults,  Na: 

Let 

Np  =  the  number  of  adult  females  resighted  or  newly  banded  prior  to 
hatching  of  eggs; 

Pb  =  the  number  of  Np  that  were  later  seen  with  broods;  and 
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B  =  the  total  number  of  marked  adult  brood  hens  sighted  in  a  given 
year. 

Assuming  that  the  same  proportion  of  brood  hens  were  among  hens  located 
prior  to  broods  hatching  (Np)  as  among  those  not  seen  prior  to  hatching 
(Na-Np)  and  that  the  same  proportion  of  brood  hens  was  found  among  Np  as 
among  Na-Np,  then 

B  =  (Pb/Np)  •  Na; 


therefore, 

Na  =  B- (Np/Pb) . 

An  analogous  calculation  yielded  an  estimated  number  of  yearlings,  N  . 

Biologically,  this  formula  will  overestimate  if  the  ratio  of  Pb/Np  is 
abnormally  low,  as  might  happen  for  example,  if,  by  capturing  a  bird,  one 
reduces  its  chances  of  successfully  nesting.  This  may  have  been  the  case 
with  some  new  bandings  in  May,  especially  yearlings  nesting  for  the  first 
time,  but  there  is  no  way,  at  present,  to  test  the  bias.  Regardless  of 
biological  biases,  this  formula  overestimates  the  number  of  females  (Ricker 
1958)  and  a  correction  could  be  added  to  compensate  for  this  overestimate. 
Corrections  were  calculated  but  they  made  little  difference  to  the  estimates. 

I  feel  that  corrections  were  unwarranted,  given  the  nature  of  the  data. 

Thus,  both  statistically  and  biologically,  this  formula  will  likely  over¬ 
estimate  the  number  of  females  and  the  estimates  of  densities  based  on  these 
calculations  are  maximums. 

I  first  tested  the  null  hypothesis  that  there  were  no  differences  in 
the  ratio  Pb/Np  between  years  or  age-classes  (Table  3).  Among  adults,  a 
greater  proportion  may  have  had  broods  in  1971  than  in  1969  (0.1>p>0.05)  but 
among  yearlings,  there  were  no  differences  between  years.  Between  age- classes, 
adults  were  more  successful  in  1971  than  yearlings  (p<0.05)  but  all  other 
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Table  3.  Data  used  to  compare  the  ratio  Pp/Np  among  adult  and  yearling 
female  blue  grouse,  t 


Year 

Adults 

Yearlings 

Number 

Proportion 

Number 

Proportion 

1968 

2/3 

0.67 

4/9 

0.45 

1969 

16/33 

0 . 49 

19/37 

0.51 

1970 

19/34 

0.56 

13/28 

0.46 

1971 

24/32 

0.75 

36/69 

0.52 

Totals 

45/69 

0.65 

72/143 

0.50 

(without  1969) 

+Pi  is  the  number  of  hens  located  before  broods  hatch  and  then  later  seen 
with  a  brood.  N  is  the  total  number  of  hens  located  before  broods  hatch. 

The  ratio  among  Sdults  in  1969  was  nearly  significantly  different  from  adults 
in  other  years  (p^  0.1),  hence  it  was  not  used  to  calculate  the  total  ratio. 
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comparisons  showed  no  differences.  Thus,  the  ratio  P^/Np>  was  estimated  as 
0.5  for  adults  in  1969  and  0.65  in  other  years.  For  yearlings  the  estimate 
was  0.5  in  all  years. 

Collection  of  Blood 

For  genetic  studies  a  sample  of  blood  was  collected  from  all  birds 
handled.  Blood  was  collected  by  clearing  about  3  cm  of  feathers  from  the 
distal  end  of  the  humerus  on  the  ventral  surface  of  the  wing  and  puncturing 
a  vein  that  crosses  the  humerus  at  this  point.  As  blood  formed  a  pool  it 
was  drawn  into  a  syringe  which  had  been  rinsed  with  4%  sodium  citrate.  A 
1  cc  disposable  tuberculin  syringe  with  27  gauge,  1/2  inch  needle  was  used 
with  small  chicks.  Otherwise,  a  2-1/2  or  3  cc  syringe  with  22  gauge  1  inch 
needle  was  used.  (Syringes  found  to  work  best  in  the  field  were  those  with 
a  heavy  plastic  case  which  prevented  damage  to  the  syringe  and  needle.)  The 
amount  of  blood  collected  varied  from  0.05  to  2  cc  depending  on  the  size  of 
the  bird.  Samples  were  successfully  collected  from  chicks  that  had  just 
hatched  as  well  as  all  older  ages.  Since  1968,  no  more  than  1  cc  of  blood 
was  taken  from  any  one  individual. 

While  in  the  field,  blood  was  stored  in  1  cc  centrifuge  tubes  on  ice 
in  a  1/2  liter  vacuum  bottle.  Upon  arrival  at  camp,  usually  within  10  hours, 
blood  was  centrifuged  and  the  plasma  was  separated  and  frozen.  Plasma  was 
usually  frozen  for  2  to  9  months  at  -20  to  -60  C.  Longer  storage  appeared 
to  have  had  little  effect  on  the  protein  examined.  I  obtained  the  same 
qualitative  results  with  14  samples  run  after  4  months'  storage  and  again 
3  years  later. 

Electrophoretic  Techniques 

All  samples  were  analyzed  for  variation  at  the  Ng  locus  (see  The  Genetic 
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Marker)  using  starch-gel  electrophoresis.  The  buffer  systems  used  for  detec¬ 
tion  of  these  bands  were  modified  from  Fujino  and  Kang  (1968).  Starch-gel 
electrophoresis  (Smithies  1955,  1959)  of  plasma  was  carried  out  in  a  discon¬ 
tinuous  system  of  buffers  (Poulik  1957).  Bridge  buffer  contained  23.6  g 
boric  acid,  4.0  g  lithium  hydroxide,  and  1  liter  of  water  (pH  7.8);  gel 
buffer  contained  2  g  citric  acid,  7  g  Tris,  170  ml  bridge  buffer,  and  1  liter 
water  (pH  8.1).  Vertical  electrophoresis  was  conducted  at  2  C  with  a  constant 
400  volts  and  a  starting  current  of  2  ma/cm  of  gel  for  about  4  hours  but  it 
was  allowed  to  proceed  until  the  borate  buffer  had  migrated  16  cm.  This 
was  important  since  the  Ng  system  was  not  clear  until  the  borate  buffer  had 
migrated  through  it. 

Gels  were  stained  for  2  hours  in  the  dark  in  a  solution  commonly  used 
for  detecting  esterases:  0.05  g  a-napthyl  acetate  (dissolved  in  2.5  ml 
acetone),  10  ml  0.1  molar  Tris-HCl  buffer  (pH  7),  and  250  mg  Fast-blue  RR 
salt  were  added  to  90  ml  water  and  mixed  for  about  two  minutes  immediately 
prior  to  staining.  Destaining  was  for  24  to  48  hours  in  methanol,  water  and 
acetic  acid  (5:5:1).  Bands  at  this  locus  were  visible  anodal  to  the  major 
albumin  esterase  activity  as  white  areas  on  a  grey  or  blue  background 
(Figure  2).  This  method  is  modified  from  that  reported  in  Birdsall  et  at . 
(1970).  All  questionable  readings  were  rerun  until  definitive  results  were 
obtained.  Less  than  1%  of  all  samples  could  not  be  classified,  but  this  was 
due  to  smearing,  not  to  an  absence  of  the  white  bands. 

Statistical  Analysis 

The  specific  statistical  tests  made  on  the  data  are  mentioned  in  the 
text.  However,  all  x2  tests  were  done  with  a  G-test  (Sokal  and  Rohlf  1969). 
Where  appropriate,  Yate's  correction  for  continuity  was  applied. 
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All  statistical  analyses  were  done  with  an  IBM/ 360  computer  using  the 
APL  language.  Most  of  the  programs  are  documented  in  Library  160  in  Armitage 
(1971).  The  level  of  significance  is  5%,  but  I  may  refer  to  trends  in  the 
data  when  the  probability  value  is  less  than  10%. 
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RESULTS 

Demography 

Demography  is  the  study  of  populations,  their  density,  age  structure, 
dispersion,  birth  and  death  rates,  and  dispersal.  Within  the  limitations 
of  this  study,  only  density  and  birth  rates  (i.e.,  breeding  success)  with 
respect  to  colonization  were  analyzed. 

Densi ty 

General  Considerations,  Before  looking  at  changes  in  density  associated 
with  colonization,  the  general  changes  in  density  are  given.  These  patterns 
may  help  interpretation  of  colonization  and  are  also  valuable  comparisons 
with  other  studies.  Since  I  consider  the  density  of  yearling  males  to  be 
the  same  as  for  yearling  females  (see  METHODS)  I  will  discuss  the  results 
for  adult  males,  adult  females  and  yearling  females  (Figure  5).  Following 
this  section  changes  in  density  associated  with  colonization  are  examined. 

Area  108e  had  a  consistently  higher  density  of  adult  males  and  females 
than  either  of  the  other  areas,  at  least  since  1969.  Area  107  had  a  higher 
density  of  adult  males  than  104,  except  in  1970,  when  the  relationship  was 
reversed.  Both  104  and  107  had  nearly  the  same  density  of  adult  females 
in  all  years. 

Recruitment  of  yearlings  in  1969  was  higher  in  108e  than  in  the  other 
areas.  In  1970  and  1971  108e  had  recruitment  rates  essentially  identical 
to  those  in  107.  In  1968,  104  had  a  higher  rate  of  recruitment  of  yearlings 
than  107,  while  in  1969  and  1970  these  two  areas  had  nearly  identical  rates 
of  recruitment.  In  1971,  107  had  a  higher  rate  of  recruitment  than  104. 

Comparisons  within  areas  but  between  years  gave  the  following  results: 

In  104  all  segments  of  the  breeding  population  declined  between  1968  and 
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Figure  5. 


Density  of  adult  male,  adult  female 
blue  grouse  on  the  three  main  study 


and  yearling  female 
areas,  1968  to  1971. 
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1969  but  after  1969  there  was  a  general  increase.  In  107,  adult  males 
increased  from  1968  to  1969,  decreased  slightly  in  1970  and  increased  markedly 
in  1971.  Adult  females  in  107  decreased  between  1968  and  1969,  and  increased 
in  both  1970  and  1971.  Recruitment  of  yearlings  into  107  increased  steadily 
since  1968  with  a  sharp  increase  in  1971.  Density  of  grouse  in  108e  was 
higher  than  in  107  but  densities  of  both  adult  males  and  females  paralleled 
those  in  107.  In  1970  and  1971,  recruitment  of  yearlings  in  108e  was  the 
same  as  in  107. 

In  general,  populations  increased  since  1969  with  increases  in  1971 
being  most  marked.  Recruitment  of  yearlings  increased  most  markedly,  especi¬ 
ally  in  108e  and  107  in  1971.  Because  of  the  high  recruitment  rates  in  1971, 
populations  in  the  Ash  River  Valley  appear  to  still  be  increasing. 

Using  these  estimates  of  density,  I  calculated  the  minimum  percentage 
of  birds  banded  in  each  area  for  each  year  (Table  4).  In  most  years  and 
areas,  greater  than  50%  of  all  sex-age  classes,  except  yearling  males,  were 
marked.  In  all  years  higher  percentages  of  adult  males  were  marked  than 
yearling  males.  This  was  due  to  behavioural  differences  between  the  two  age 
groups.  The  percentage  of  females  banded  was  as  high  or  higher  than  for 
adult  males,  suggesting  I  was  equally  successful  at  locating  and  marking 
females  and  adult  males.  The  low  percentage  of  adult  males  tagged  in  104 
in  1971  (42%)  was  probably  due  to  two  factors:  1)  The  increase  in  total 
density  of  birds  on  all  study  areas  in  1971  made  complete  coverage  of  all 
areas  difficult,  especially  early  in  spring  and  early  in  the  year  less 
effort  was  expended  on  104  than  on  107  or  108e.  2)  Many  of  the  birds  on  104 

seemed  wild  in  1971,  making  resighting  and  capture  difficult.  The  wildness 
was  found  in  birds  on  104  in  other  years  but  it  was  more  distinct  in  1971. 

Changes  in  Density  Associated  with  Colonization.  The  data  on  density 


. 
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Table  4.  Percentage  of  blue  grouse  banded  on  the  study  areas,  1968  to  1971. 

The  percentage  of  adult  males  is  based  on  the  number  of  territories 
counted  and  the  number  of  marked  territory  holders.  The  yearling 
male  percentage  was  based  on  the  estimated  number  of  females  (year¬ 
ling)  and  the  actual  number  of  yearling  males  banded.  Female  per¬ 
centages  were  based  on  the  actual  number  of  females  marked  and  the 
total  number  estimated  to  be  alive.  The  numbers  in  brackets  represent 
actual  number  banded/total  estimated  number. 


Area  Age  Sex  1968  1969  1970  1971 


104 

Adult 

4 

46 

(13/28) 

70 

(23/33) 

75 

(28/37) 

42 

(21/50) 

Yearling 

15 

(  5/34) 

46 

(15/32) 

17 

(  8/48) 

30 

(17/58) 

Adult 

? 

72 

(23/32) 

64 

(27/42) 

70 

(32/46) 

71 

(34/48) 

Yearling 

? 

71 

(24/34) 

85 

(27/32) 

54 

(26/48) 

57 

(33/58) 

107 

Adult 

55 

(17/31) 

59 

(25/42) 

70 

(32/46) 

70 

(44/63) 

Yearling 

d 

28 

(  5/18) 

50 

(15/30) 

29 

(14/49) 

38 

(42/101) 

Adult 

? 

66 

(24/36) 

88 

(30/34) 

77 

(44/57) 

72 

(41/57) 

Yearling 

? 

78 

(14/18) 

100 

(30/30) 

71 

(35/49) 

76 

(81/101) 

108e 

Adult 

d1 

— 

86 

(19/22) 

90 

(17/19) 

82 

(22/27) 

Yearling 

-- 

31 

(  9/29) 

29 

(  4/14) 

29 

(10/35) 

Adult 

¥ 

-- 

68 

(22/32) 

94 

(30/32) 

76 

(29/38) 

Yearling 

¥ 

-- 

80 

(23/29) 

72 

(10/14) 

69 

(24/35) 

. 
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presented  in  the  previous  section  did  not  show  changes  occurring  on  differ¬ 
ent  age  habitats.  Since  a  mosaic  of  habitat  types  existed  within  each  study 
area,  and  since  each  habitat  type  represented  different  years  since  logging, 
it  is  of  interest  to  know  if  birds  were  distributed  at  random  over  these 
various  habitat  types,  or  whether  they  were  more  dense  on  some  areas  than 
others . 

In  order  to  calculate  density  of  birds  on  each  age  habitat,  I  calculated 
the  size  of  each  habitat  (Table  5),  the  number  of  resightings  on  that  habitat 
(Table  6),  and  the  percentage  of  the  total  number  of  birds  banded  (Table  4). 
Samples  for  each  area  (104,  107  and  108e)  were  too  small  for  between-area 
comparisons,  hence  data  for  all  areas  were  combined.  Even  with  data  com¬ 
bined  between  areas,  samples  were  usually  too  variable  for  between-year 
comparisons  and  all  years  were  combined  for  final  analysis  (Table  7). 

Density  of  adult  males  and  adult  and  yearling  females  for  all  age- 
habitats  are  given  in  Figure  6.  Density  of  adult  males  and  females  increased 
steadily  from  one  through  at  least  seven  years  after  logging.  Age-class  >8 
is  composed  of  all  habitat  8  years  old  and  older  and  includes  some  habitats 
which  were  unsuitable  (according  to  Bendell  and  Elliott  1966).  Regression 
on  these  points  (excluding  point  £8)  showed  that  adult  males  and  adult  females 
were  increasing  at  about  2.9  to  3.4  birds/km2/age  of  habitat  during  this 
seven  year  period  (Yma^es  =  l* 52+2. 89X  [p^O.Ol]  ;  ^females  =  1*71+3. 38X 
[psO.Ol])  while  density  of  yearling  females  did  not  increase  as  rapidly  as 
adults  (Y  =  7*51+2*23X  [p^ 0.01]) .  These  density  figures  show  a  sharp 
increase  in  density  of  adults  and  smaller  increases  in  density  of  yearlings 
in  the  first  seven  years  following  logging. 

Density  of  yearling  females  did  not  fit  a  linear  model  well.  There 
appeared  to  be,  among  yearlings,  a  two  stage  period  of  colonization. 
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Table  5.  Size  (in  km2)  of  each  age  habitat  studied,  1968  to  1971,  for 
all  areas  combined. 


Year 

Age  of  ' 

Habitat 

1 

2 

3 

4 

5 

6 

7 

>8 

1968 

0.63 

0.82 

1.07 

0.72 

0.41 

0.12 

1.77 

1969 

1.20 

0.76 

1.40 

1.10 

1.20 

-- 

0.41 

1.96 

1970 

0,82 

1.20 

0.76 

1.40 

1.10 

1.20 

— 

2.26 

1971 

0.27 

0.82 

1.20 

0.76 

1.40 

1.10 

1.20 

2.26 

Totals 

2.92 

3.60 

4.43 

3.98 

3.70 

2.71 

1.73 

8.25 

■ 
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Table  6. 

Number  of  blue  grouse  captured  or 
habitat,  1968  to  1971. 

regi 

stered  on 

each 

age  of 

Year 

Age 

Sex 

Age 

:  of  Habitat 

Total 

1 

2 

3 

4 

5 

6 

7 

>8 

1968 

Adult 

d" 

0 

3 

5 

6 

6 

0 

10 

30 

Yearling 

<? 

0 

0 

2 

3 

2 

0 

3 

10 

Adult 

$ 

0 

3 

5 

11 

-- 

4 

0 

24 

47 

Yearling 

o 

i - 

1 

6 

6 

7 

__ 

10 

0 

8 

58 

1969 

Adult 

4 

6 

14 

8 

18 

— 

7 

10 

67 

Yearling 

d1 

5 

3 

6 

2 

13 

— 

7 

4 

40 

Adult 

? 

/L 

r\ 

6 

14 

11 

24 

-r- 

8 

12 

79 

Yearling 

9 

4 

12 

15 

9 

21 

4 

16 

81 

1970 

Adul  t 

d" 

1 

7 

4 

19 

7 

17 

— 

22 

77 

Yearling 

d 

3 

5 

3 

5 

0 

4 

— 

6 

26 

Adult 

? 

2 

9 

13 

21 

9 

24 

-- 

28 

106 

Yearling 

? 

5 

7 

6 

12 

10 

13 

-- 

18 

71 

1971 

Adult 

d 

6 

0 

4 

4 

17 

9 

22 

25 

87 

Yearling 

d 

9 

6 

4 

8 

8 

3 

11 

19 

68 

Adult 

? 

0 

9 

8 

9 

19 

9 

26 

25 

105 

Yearling 

? 

11 

16 

10 

9 

30 

17 

24 

22 

139 

Total 

Adult 

d 

11 

16 

27 

37 

42 

32 

29 

67 

261 

Yearling 

d 

17 

14 

15 

18 

21 

9 

18 

32 

144 

Adult 

? 

4 

29 

40 

52 

52 

37 

34 

89 

337 

Yearling 

? 

21 

41 

37 

37 

61 

40 

28 

64 

329 

30 


Table  7.  Density  of  blue  grouse/km2  on  various  age  habitats  for  1968  to 
1971  combined.  (Yearling  males  are  considered  the  same  as 
yearling  females.) 


Years 


Age  of  Habitat 


1 

2 

3 

4 

5 

6 

7 

£8 

Adult 

d1 

5.7 

6.2 

9.7 

12.7 

16.6 

18.2 

22.5 

13.9 

Adult 

? 

1.8 

10.8 

14.3 

15.6 

19.7 

17.3 

27.3 

14.8 

Yearling 

? 

9.8 

13.9 

12.9 

12.2 

22.4 

21.7 

22.2 

11.4 
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Figure  6.  Density  of  adult  male,  adult  female  and  yearling  female 
blue  grouse  on  each  age  habitat  for  all  years  combined. 
Age  of  habitat  is  equivalent  to  years  since  logging. 
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Populations  were  fairly  stable  at  about  10-13  yearling  females/km2  during 
the  first  four  years  after  logging  but  five  years  after  logging  populations 
increased  about  70%  to  a  new  level  of  about  22  yearling  females/km2.  This 
new  level  was  stable  for  at  least  three  years. 

All  groups  declined  in  age-class  <s8.  Since  age-class  ^8  was  a  composite 
and  included  some  habitat  which,  according  to  Bendell  and  Elliott  (1966), 
could  be  called  unsuitable,  it  is  difficult  to  say  what  happened  in  older, 
but  suitable,  habitat.  However,  densities  probably  stabilized  shortly  after 
year  seven. 

Age-specific  density  and  age  structure  of  grouse  living  on  each  habitat 
type  were  closely  related.  For  analysis  of  age  structure,  only  birds  known 
to  be  alive  were  used.  As  the  data  on  density  would  suggest,  there  was  a 
high  proportion  of  yearlings  on  age-classes  1  and  2  (Figure  7).  Age  structure 
of  birds  on  one-year-old  habitat  was  about  70%  yearlings  and  was  significantly 
different  from  age  structure  of  birds  on  all  other  habitat  types  (p$10~5). 

Age  structure  of  birds  on  two-year-old  habitat  was  about  55%  yearlings  and 
was  significantly  different  from  age  structure  of  birds  on  older  habitat 
types  (p^O.Ol).  Age  structure  relationships  on  all  other  habitat  types  show 
no  differences  or  trends,  suggesting  age  structure  stabilized  at  about  45% 
yearlings  by  three  years  after  logging. 

Recruitment  onto  one-  and  two-year- old  habitat  was  related  with  produc¬ 
tion  the  previous  summer.  For  example,  1970  was  a  year  of  very  high  repro¬ 
ductive  output  and  in  1971  recruitment  onto  one-year-old  habitat  was  higher 
than  had  been  recorded  previously.  The  total  number  of  yearlings  captured 
in  1971  was  nearly  double  the  number  captured  in  any  previous  year,  but  x2 
analysis  suggested  that  recruitment  was  significantly  higher  in  one-  and  two- 
year  old  habitat  than  in  other  habitat  types.  In  order  to  test  the  null 
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Figure  7.  Percentage  of  yearling  male  and  female  blue  grouse  among 

all  banded  grouse  on  each  age  habitat,  all  years  combined. 
The  male  percentage  is  lower  than  the  female  percentage 
because  yearling  males  are  more  difficult  to  locate  than 
yearling  females. 


PERCENT  YEARLINGS 


AGE  OF  HABITAT 
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hypothesis  of  no  differences  in  recruitment  into  all  age  habitats  between 
years,  x2  goodness-of-fit  values  were  calculated.  These  calculations  show 
number  of  yearlings  on  age-class  1  and  2  habitat  in  1971  was  significantly 
larger  than  in  the  other  three  years  (Table  8).  Other  habitat  types  had 
smaller  values,  except  age-class  *8. 

In  summary,  density  of  adult  males  and  females  increased  for  at  least 
seven  years  following  logging.  Yearlings  occurred  in  higher  absolute 
densities  than  adults  on  younger  age  habitats  but  did  not  increase  as  rapidly 
as  adults.  Yearlings  seemed  to  have  a  two  stage  pattern  of  settlement.  The 
first  stage  was  at  about  10  to  12  yearling  females/km2  and  lasted  four  years. 
The  second  stage  was  at  about  22  yearling  females/km2  and  lasted  at  least 
three  years.  All  segments  of  the  population  declined  in  habitat  ^8  years  old 
but  this  may  have  been  due  to  some  of  this  old  habitat  being  unsuitable. 

There  were  signifcant  shifts  in  age  structure  of  colonizing  populations 
from  about  70%  to  45%  yearlings  in  three  years.  Recruitment  of  yearlings 
onto  young  habitat  (one  and  two  years  old)  was  significantly  higher  in  1971, 
a  year  following  very  high  production,  than  in  other  years. 

Breeding  Success 

General  Considerations .  Several  calculations  associated  with  breeding 
may  be  needed  for  an  evaluation  of  population  changes.  Three  of  these  are 
the  percentage  of  females  with  broods,  the  average  size  of  broods  and  the 
total  estimated  chick  production  (which  may  be  important  in  evaluating  the 
number  of  colonizers  available  for  the  following  year). 

I  assume  all  females  were  capable  of  and  attempted  to  breed  each  year 
(Zwickel  and  Bendell  1967).  However,  not  all  females  were  successful  at 
completing  the  nesting  cycle.  The  proportion  which  was  successful  was  part 
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Table  8.  Results  of  x2  analyses  of  captures  of  yearling  blue  grouse  on 
various  age  habitats  in  1971.  A  significant  departure  means 
that  there  were  a  greater  number  captured  in  that  age  habitat 
in  1971  when  compared  to  the  years  1968  to  1970  combined. 


Age  of 
Habitat 


Yearling  Males  Yearling  Females 


Totals 


39.30 

0.00 

47.3 

0.00 

86.50 

0.00 

3.18 

0.07 

6.11 

0.01 

9.20 

0.00 

0.01 

0.90 

0.01 

0.92 

0.02 

0.88 

7.34 

0.01 

0.61 

0.44 

4.79 

0.02 

0.00 

0.98 

3.02 

0.08 

2.21 

0.14 

0.18 

0.67 

0.09 

0.77 

0.01 

0.93 

0.17 

0.68 

4.96 

0.03 

2.19 

0.14 

15.00 

0.00 

1.53 

0.22 

10.70 

0.00 

8 
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of  the  effective  breeding  number  (Kimura  and  Crow  1963).  Effective  breeding 
was  analyzed  by  two  methods:  First,  I  calculated  the  total  percentage  of 
all  females  with  broods.  This  overestimated  the  percentage  with  broods 
since  successful  hens  with  broods  were  easier  to  locate  than  unsuccessful, 
lone  hens.  About  50%  of  the  hens  were  first  located  after  broods  hatched. 
Thus,  this  calculation  likely  gave  maximum  estimates.  Alternatively,  the 
proportion  of  females  found  prior  to  broods  hatching  that  later  had  a  brood 
was  calculated.  This  estimate  was  a  minimum  and  was  previously  calculated 
in  connection  with  estimating  numbers  (Table  3).  These  two  estimates  set 
limits  on  what  I  called  effective  breeding. 

Both  sets  of  data  are  presented  in  Table  9.  Data  for  1968  are  weak 
since  at  that  time  I  did  not  spend  the  necessary  time  on  the  intensively 
worked  study  areas,  especially  early  in  the  year,  to  make  large  enough 
counts  of  females.  Thus,  most  of  my  conclusions  are  based  on  data  for  the 
latter  three  years. 

In  most  years  between  50  and  90%  of  adults  and  between  45  and  85%  of 
yearlings  were  effective  breeders.  There  seemed  to  be  a  smaller  percentage 
of  effective  breeding  among  adults  in  1969  than  in  either  1970  or  1971. 
Between  age-classes  and  within  years,  adults  were  more  successful  than  year¬ 
lings  only  in  1971. 

Zwickel  and  Bendell  (1967a)  suggested  a  maximum  of  50%  of  females  were 
broodless  in  a  given  year  in  their  study.  Their  calculations  were  based  on 
total  counts  of  females  on  breeding  range.  Using  these  criteria  (maximum 
estimates),  my  data  suggest  that  less  than  30%  of  females  were  unsuccessful 
in  a  given  year.  The  difference  between  the  two  studies  was  significant  and 
two  explanations  are  possible:  First,  perhaps  I  was  more  efficient  at 
finding  females,  especially  females  with  broods.  I  used  an  imitation  chick 
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Table  9.  Maximum  and  minimum  percentages  of  female  blue  grouse  with 

broods  on  all  study  areas  combined,  1968  to  1971.  The  minimum 
percentages  are  from  Table  3  and  the  maximum  percentages  are 
based  on  the  method  of  Zwickel  and  Bendell  (1967a). 


Year 

Age 

Actual  No. 
of  Females 

Actual  No. 
of  Broods 

%  with 

Maximum 

Broods 

Minimum 

1968 

Adult 

46 

44 

96 

67 

Yearling 

38 

26 

69 

45 

1969 

Adult 

79 

55 

68 

49 

Yearling 

80 

46 

56 

51 

1970 

Adult 

106 

88 

83 

56 

Yearling 

67 

53 

79 

46 

1971 

Adult 

104 

95 

91 

75 

Yearling 

138 

100 

73 

52 
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call  extensively  during  searches  once  broods  hatched  and  this  caused  many 
females  to  respond,  thus  identifying  their  location.  Zwickel  (pers.  oomm. ) 
did  not  use  this  call.  Alternatively,  perhaps  this  was  a  case  of  selection 
for  genotypes  with  a  high  reproductive  rate  in  an  expanding  population.  In 
the  study  of  Zwickel  and  Bendell  (1967a),  the  population  was  relatively 
stable,  whereas  populations  which  I  studied  increased  rapidly  after  1969. 
Selection  for  high  and  early  reproduction  in  rapidly  expanding  populations 
has  been  hypothesized  on  theoretical  grounds  (Cole  1954,  Lewontin  1965, 
MacArthur  and  Wilson  1967). 

Among  females  with  broods,  I  calculated  average  brood  size  by  counting 
the  number  of  chicks  in  broods  known  to  be  older  than  15  days  [by  backdating 
of  captured  chicks  (Zwickel  and  Lance  1966)]  or  estimated  older  than  20  days 
(by  size  of  chicks).  A  one-way  analysis  of  variance  with  unequal  sample 
sizes  was  used  to  ascertain  if  there  were  differences  between  parental  age- 
classes,  years,  or  areas  (Steel  and  Torrie  1960,  Sokal  and  Rohlf  1969). 

The  means  and  standard  errors  for  all  counts  are  presented  in  Table  10. 

Adults  had  significantly  larger  brood  sizes  than  yearlings  in  all  years 
except  1968  (1968:  F=0.04,  df=l,66,  p^0.75;  1969:  F=5.23,  df=l,95,  p^0.05; 

1969:  F=5 . 23,  df=l,9S,  ps0.05;  1970:  F=16.56,  df=l,164,  psO.OOl;  1971: 

F=13.14,  df=l,164,  p^O.OOl).  Thus,  parental  ages  were  considered  separately 
in  further  analyses.  Brood  sizes  of  adults  between  areas  showed  no  signifi¬ 
cant  differences  (1968:  F=0.03,  df=l,41,  p$0.75;  1969:  F=0.22,  df=2,50, 

psO.75;  1970:  F=2.63,  df=2,99,  p*0.1;  1971:  F=0.67,  df=2,90,  p$0.75). 

Brood  sizes  of  yearlings  between  areas  showed  significant  differences  only 
in  1968  (1968:  F=6.89,  df=l,23,  p*0.001;  1969:  F=1.09,  df=2,41,  p;>0.5; 

1970:  F=0. 25 ,  df=2.61,  ps0.75;  1971:  F=0.93,  df=2,70,  psO.5).  In  1968  the 

difference  stemmed  from  a  homogeneous  set  of  six  brood  counts  and  I  disregarded 
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Table  10.  Average  (±1  standard  error)  brood  sizes  for  adult  and  yearling 
female  blue  grouse,  1968  to  1971,  in  all  study  areas. 


Year 


Age 


104 


Area 


107 


lOSe 


Totals 


1968 


Adult 


N 


21 

3 . 5  ±0 . 4 


22 

3 . 6  ±0 . 3 


43 

3.510.2 


Yearling  N 


19 

4.2±0.7 


6 

1.711.2 


9  ^ 


3.610.8 


1969 


Adult 


N  20 

3.510.5 


2'0 

3.810.5 


13 

3.911.2 


55 

3.710.4 


Yearling  N 


11 

2.911 .0 


23 

3.210.3 


10 

2.411 .0 


41 

2.910.4 


1970 


Adult 


N 


33 

4.210.6 


44 

5.010.2 


25 

5.211.0 


102 

4.810.3 


Yearling  N  29 


3.810.6 


30 

3.710.6 


3.214.6 


64 

3.710.5 


1971 


Adult 


N  24 

5.H1.0 


40 

4.710.4 


29 

4.410.6 


93 

4.710.4 


Yearling  N  23  36  14  73 

4.010.8  3.310.3  3.611.5  3.610.4 


1 
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this  difference.  Areas  were  combined  and  between-year  analyses  calculated. 
These  calculations  showed  highly  significant  between-year  fluctuations  in 
average  brood  size  of  adults  (F=8.04,  df=3,287,  p^O.OOl)  but  not  for  year¬ 
lings  (F=l . 96 ,  df=3 , 202 ,  psO.l). 

The  between-year  difference  among  adults  was  reduced  to  a  difference 
between  1968+1969  and  1970+1971.  Adults  in  these  latter  two  years  had 
large  brood  sizes  for  grouse  on  Vancouver  Island.  Thus,  average  brood  size 
of  adults  seemed  more  variable  than  for  yearlings. 

Average  brood  sizes  for  both  yearlings  (2.9  to  3.7)  and  adults  (3.5 
to  4.8)  was  high  when  compared  to  data  from  other  studies  on  Vancouver 
Island.  Bendell  (1955)  reported  an  average  brood  size  of  2  to  2.7  in  late 
summer,  1950  to  1952.  Zwickel  and  Bendell  (1967a)  reported  an  average  brood 
size  of  2.2  to  2.9  for  late  summer  during  1962  to  1964  (adults  and  yearlings 
combined).  It  is  not  easy  to  explain  this  difference  as  artifact.  Once 
chicks  can  fly  (at  about  7  to  10  days  of  age),  they  were  relatively  easy  to 
count  and  both  Zwickel  (pers.  comm.)  and  I  feel  that  our  counts  were  accur¬ 
ate.  A  possible  explanation  is  that,  as  for  the  percentages  of  females 
breeding,  there  was  recent  selection  for  high  reproductive  genotypes  in  the 
population.  Whether  this  hypothesized  selection  was  the  cause  of  population 
increases  or  the  result  of  them  is  not  known.  In  support  of  this  hypothesis, 
Zwickel  reopened  his  studies  of  blue  grouse  on  Vancouver  Island  and  found 
that  the  population  at  his  old  study  area  was  increasing;  the  average  brood 
size  was  also  higher  than  in  his  earlier  studies  (Zwickel,  pers.  eomm.). 

Zwickel  and  Bendell  (1967a)  concluded  that  recruitment  in  a  stable 
population  was  not  dependent  on  the  previous  year’s  production.  I  compared 
these  two  variables  by  estimating  total  production  on  each  study  area  in  a 
given  year  and  calculating  correlation  coefficients  between  production  and 
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number  of  yearlings  captured  in  the  following  year  (Table  11).  There  was 
a  nearly  significant  correlation  between  these  two  variables  for  all  data 
combined  (r=0.970,  p<0.2)  and  for  areas  considered  separately  (r=0.679, 
p-0.05)  suggesting  that  recruitment  in  an  expanding  population  may  have 
been  related  to  the  previous  year's  production. 

Summarizing,  a  high  percentage  of  females  was  successfully  breeding 
in  the  Ash  River  Valley.  Adults  were  slightly  more  successful  than  yearlings. 
Among  successful  hens,  average  brood  size  was  also  high.  Adults  had  larger 
broods  than  yearlings  in  all  years  except  1968  but  there  was  little  varia¬ 
tion  between  areas.  Brood  sizes  of  adults,  but  not  yearlings,  varied 
significantly  from  year  to  year.  The  variation  was  reduced  to  a  difference 
between  1968+1969  and  1970+1971.  In  the  latter  two  years  brood  sizes  of 
adults  were  large.  Correlation  between  recruitment  in  one  year  and  produc¬ 
tion  in  the  previous  one  was  found.  The  data  support  the  hypothesis  that 
in  expanding  populations  there  is  selection  for  high  reproductive  rate. 

Breeding  Success  and  Colonization.  In  the  previous  section  I  presented 
evidence  suggesting  higher  breeding  success  in  expanding  populations  of 
grouse.  Now  I  ask  the  question:  Does  breeding  success  vary  with  birds  on 
different  age  habitats?  Since  population  densities  were  increasing  with 
age  of  habitat  it  was  reasonable  to  suppose  that  there  would  be  a  systematic 
change  in  breeding  success.  As  in  the  previous  section,  percentage  of 
females  with  chicks  and  average  brood  size  were  used  as  a  measure  of  breeding 
success . 

Since  data  were  too  few  to  compare  minimum  estimates,  I  compared  maxi¬ 
mum  percentages  of  females  known  to  have  had  broods  for  adults  and  yearlings 
on  different  age-habitats  (Table  12).  Females  on  older  age-habitats  tended 
to  have  a  higher  percentage  with  broods  than  females  on  younger  age  habitats. 
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This  adjustment  is  based  on  a  predicted  percentage  contribution  from  area  108e  based  on  the  contributions 
in  the  following  two  years. 
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Table  12.  Maximum  proportion  of  female  blue  grouse  on  each  age  habitat  with  and  without  broods. 
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I  first  tested  the  null,  hypothesis  that  there  were  no  differences  between 
adults  and  yearlings  within  a  habitat,  type.  Only  age-class  7  had  a  signi¬ 
ficantly  higher  proportion  of  adults  than  yearlings  with  broods  (p£0.03). 

Since  there  were  no  other  obvious  trends  in  these  data,  all  data  were 
combined.  When  x2  was  calculated  on  pairs  of  samples,  the  results  indicated 
that  age-class  7  and  b8  had  a  significantly  larger  percentage  of  successfully 
breeding  females  than  other  age  habitats. 

These  results  are  potentially  biased,  however.  If  in  older  habitats 
it  was  more  difficult  to  resight  a  female  before  broods  hatch  due  to  the 
structure  of  vegetation,  then  the  estimated  percentage  of  females  with 
broods  would  be  too  high.  Secondly,  Mos.sop  (1971)  suggested  that  birds  in 
populations  at  high  densities  were  wilder  than  birds  in  populations  at  low 
densities.  If  true,  then  it  may  have  been  more  difficult  to  resight  females 
without  chicks,  and  this  would  increase  the  observed  percentages  with  broods. 
At  any  rate,  the  data  do  not  support  the  hypothesis  that  birds  in  low  density 
populations  would  be  more  successful  breeders  than  birds  in  higher  density 
populations . 

I  calculated  the  average  number  of  chicks  per  brood  hen  on  each  age 
habitat  (Table  13).  There  were  no  significant  differences  in  the  relation¬ 
ship  between  age  of  habitat  and  brood  size  in  1968,  1969,  or  1971.  In  1970, 
among  yearling  females  there  were  significant  differences  in  brood  sizes 
between  habitat  types.  Although  there  were  differences,  there  were  no  trends 
in  the  data,  with  brood  sizes  of  females  on  2-,  4-  and  8-year-old  habitat 
being  larger  than  those  on  1- ,  3-,  5-  and  6-year-old  habitat  (F=4.18,  df-6,57, 
p^0.005).  These  data  indicate  that  number  of  chicks  raised  per  hen  was 
independent  of  age  of  habitat. 

Both  the  data  for  percentage  of  females  breeding  and  number  of  chicks 


Table  13.  Average  size  of  broods  for  blue  grouse  located  on  each  age  habitat,  1969  to  1971.  Brood  coun 
include  broods  older  than  15  days  of  age.  Sample  sizes  in  1968  were  too  small  to  consider  in 
this  analysis.  F  is  the  result  of  analysis  of  variance  between  ages  of  habitats;  N  is  sample 
size;  x  is  mean;  and  se  is  the  standard  error  of  the  mean. 
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per  successful  .female  show  no  increases  in  younger  age  habitats  as  hypothe¬ 
sized.  These  results  do  not  necessarily  negate  the  hypothesis  as  stated, 
however,  since  the  population  has  been  increasing  since  1969.  There  was  a 
high  average  brood  size  and  a  high  percentage  of  females  successfully 
breeding  throughout  these  increasing  populations,  and  this  may  be  over¬ 
riding  the  hypothesis.  As  noted  earlier,  both  the  average  brood  size  and 
percentage  of  females  successful  at  ■pp^t'ing  was  larger  in  this  study  than 
reported  for  other  recent  studies  on  Vancouver  Island.  Thus,  when  these 
data  from  expanding  populations  are  compared  to  data  from  other,  more  stable, 
populations  the  general  hypothesis  was  substantiated. 

Genetics  of  the  Ng  Locus 

The  Ng  locus  was  previously  described  only  in  a  note  (Birdsall,  et  at. 
1970),  Therefore,  it  is  appropriate  to  discuss  the  inheritance  of  the  bands 
detected,  the  spatial  and  temporal  distribution  of  the  alleles  and  the  fit¬ 
ness  of  the  genotypes  before  considering  changes  in  this  system  associated 
with  colonization. 

Inheritance  of  the  bands  is  essential  since,  in  order  to  test  a  genetic 
hypothesis,  one  must  be  sure  he  is  working  with  a  genetically  determined 
character.  Inheritance  was  studied  by  comparing  the  pattern  of  variability 
to  other  systems  with  similar  characteristics  and  known  inheritance,  and 
also  by  the  method  of  incomplete  family  data  (Cooper  1968);  i.e.,  where  the 
genotype  of  one  parent  and  one  or  more  offspring  are  known. 

Understanding  spatial  and  temporal  distribution  of  the  alleles  is 
necessary,  since  potentially  the  most  useful  genetic  markers  in  the  study  of 
interaction  between  genetics  and  demography  are  those  which  are  spatially 
widespread  and  temporally  stable.  Spatial  and  temporal  variation  was  studied 
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by  comparing  frequencies  of  alleles  and  genotypes  between  samples  from 
different  parts  of  Vancouver  Island  and  from  the  adjacent  mainland,  by 
comparing  allele  frequencies  between  years  within  areas,  and  by  Hardy- 
Weinberg  calculations  [see  Falconer  (1960)  or  Li  (1955)  for  discussion  of 
implications  of  Hardy- Weinberg  calculations].  Fitness  calculations  are 
needed  to  help  understand  if  and  how  the  various  genotypes  may  be 
affecting  genetic  structure  and  demographic  events.  Also,  fitness  calcu¬ 
lations  may  shed  some  light  on  the  mechanism  of  maintenance  of  the  poly¬ 
morphism.  Reproductive  success  and  survival  of  different  genotypes  was 
used  to  calculate  fitness  of  various  genotypes. 

Inheri tance 

General  Considerations,  Three  white  bands  were  detected  after  electro¬ 
phoresis  of  plasma  from  blue  grouse.  All  individuals  scored  had  at  least 
one  band,  but  none  had  more  than  two  bands.  In  addition,  plasma  from  the 
same  birds  bled  at  two  different  times  produced  the  same  pattern.  Since 
bands  with  similar  characteristics  have  been  found  to  be  inherited  as  co¬ 
dominant  alleles  at  a  single  autosomal  locus  (Manwell  and  Baker  1970)  and 
since  there  was  close  agreement  to  Hardy- Weinberg  expectations  based  on  this 
hypothesis,  it  was  reasonable  to  suppose  that  these  bands  were  products  of 
three  codominant,  autosomal  alleles,  NgF 3  NgM,  and  NgS  for  the  fastest, 
intermediate  and  slowest  band,  respective-ly  (Figure  2,  p.7).  All  six  geno¬ 
types  from  these  three  alleles  were  found.  Homozygotes  were  characterized 
by  having  a  single  band  and  heterozygotes  by  having  two  bands. 

NgM  was  by  far  the  most  common  allele  with  a  frequency  of  near  0.8  in 
most  populations.  Because  of  the  rarity  of  NgFy  it  was  considered  equivalent 
to  NgS  for  some  analyses.  Most  populations  examined  were  in  Hardy-Weinberg 
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equilibrium  and  the  three  alleles  were  widespread  in  blue  grouse  on  Vancouver 
Island  (see  Hardy-Wei nberg  Calculations  and  Distribution).  Thus,  the  poly¬ 
morphism  appeared  stable  (Ford  1965a) . 

Incomplete  Family  Data.  Incomplete  family  data  (IFD)  were  used  to 
examine  the  genetics  of  the  Ng  locus.  IFD  are  genetic  data  where  the  geno¬ 
type  of  one  parent  and  one  or  more  offspring  are  known.  Cooper  (1966,  1968) 
and  Cooper  and  Rendel  (1968a)  pointed  out  the  utility  of  IFD  in  the  analysis 
of  genetics  of  populations.  IFD  have  the  advantage  over  population  data  in 
that  they  allow  one  to  test  hypotheses  about  inheritance  and  detect  the 
effects  of  selection  and  population  subdivision  (the  Wahlund  effect,  Wahlund 
1928).  Additionally,  IFD  are  easier  to  collect  than  complete  family  data 
and  may  be  the  only  data  available  from  most  natural  populations. 

Cooper  (1968)  urged  the  collection  of  IFD  from  marsupials  and  domestic 
species  of  mammals.  Manwell  and  Baker  (1970,  p.24)  state  that  IFD  can  only 
be  collected  from  viviparous  or  ovoviviparous  vertebrates.  These  authors 
seem  to  have  overlooked  the  possible  source  of  information  obtainable  from 
birds.  In  many  birds,  including  blue  grouse,  a  mother  and  young  are  closely 
associated  for  some  time  which  leads  to  the  possible  collection  of  IFD. 

The  method  using  IFD  is  based  on  the  fact  that  among  parents  of  known 
genotypes  there  are  certain  expectations  for  the  ratios  of  offspring.  These 
expectations  are  presented  in  Table  14,  as  an  idealized  IFD  table  with  no 
selection  or  population  subdivision.  In  this  table  there  are  two  alleles, 

A  and  B,  at  a  single  autosomal  locus,  with  frequencies  p  and  q  in  the  parent 
generation.  Offspring  are  represented  as  0  and  subscripts  are  used  to 
denote  parent  genotype  and  offspring  genotype,  respectively.  Thus,  0A  ^ 
denotes  an  AB  offspring  born  from  an  AA  parent.  From  this  table  it  can  be 
shown  that  the  following  conditions  hold: 
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Table  14.  Idealized  incomplete  family  data  table  when  there  is  no 

population  subdivision  or  selection.  0  represents  offspring. 
Subscripts  appended  to  0  signify  parental  and  offspring 
genotype  respectively.  This  table  is  for  two  allele  case 
only. 


Parental  Genotypes 
and  their  Frequencies 

Frequencies 

of  Genotypes 

among  Offspring 

Total 

Female 

Male 

AA 

AB 

BB 

AA 

AA 

P2 

P2 

p4 

- 

- 

AA 

AB 

p2 

2pq 

p3q 

p3q 

- 

AA 

BB 

P2 

q2 

- 

p2q2 

AA 

9 

p2 

p3 

p2q 

P2 

[Oa.a] 

[°A.AB] 

[°A.  ] 

AB 

AA 

2pq 

P2 

p3q 

p3q 

- 

AB 

AB 

2pq 

2pq 

p2q2 

2p2q2 

PV 

AB 

BB 

2pq 

q2 

- 

pq3 

pq3 

AB 

? 

2pq 

p2q 

pq 

pq2 

2pq 

i — i 

o 

> 

CO 

• 

> 

1 _ 1 

[°ab.ab] 

[°AB.B] 

[°AB. ] 

BB 

AA 

q2 

p2 

- 

p2q2 

- 

BB 

AB 

q2 

2pq 

- 

pq3 

pq3 

BB 

BB 

q2 

q2 

- 

- 

q4 

BB 

9 

q2 

- 

pq2 

q3 

q2 

[°B.AB] 

[°b.b] 

[°B.  ] 

Total 

p2 

2pq 

q2 

1 

[O.A] 

[O.AB] 

[O.B] 

. 
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AB-AB 


1 

1 


The  first  ratio  is  the  1:1  ratio  and  the  next  three  the  gametic  gene  ratios 


(Cooper  1968).  Using  these  ratios,  several  different  x2  goodness-of-fit 


tests  can  be  carried  out  on  data  collected  from  blue  grouse  to  detect  selec¬ 
tion  or  population  subdivision. 

The  1:1  ratio  will  only  be  affected  by  selection  operating  on  zygotes, 
while  the  gametic  gene  ratios  will  be  affected  by  gametic  selection,  zygotic 
selection,  and  population  subdivision  (Cooper  1968). 

IFD  collected  from  blue  grouse  from  1969  to  1971  are  presented  in  Table 
15.  No  other  breeding  data  in  blue  grouse  were  available.  Cooper  (1968) 
pointed  out  that  in  calculating  heterogeneity  x2  values  for  the  gametic  gene 
ratios,  the  usual  x2  test  must  be  corrected  for  multiparity  (Table  16).  The 
corrected  formula  for  heterogeneity  x2  is  given  in  Cooper  (1968) . 

Frequency  of  genotypes  among  mothers  showed  no  significant  departures 
from  Hardy-Weinberg  expectations  (Table  17).  Frequency  of  genotypes  among 
chicks,  however,  had  two  significant  departures:  in  1969  there  was  a  signi¬ 
ficant  deficiency  of  heterozygotes  and  in  1970  there  was  a  significant  excess 
(Table  17).  There  were  no  departures  from  expectations  in  1971,  nor  were 
there  any  significant  departures  from  expectations  when  the  data  were  combined 
for  all  years.  The  two  deviations  found  in  chicks  in  1969  and  1970  cancelled 
each  other  so  that  totals  for  chicks  were  in  equilibrium. 


Results  of  tests  made  on  IFD  for  agreement  with  the  1:1  ratio  are  given 
in  Table  18.  In  all  years  heterozygous  mothers  produced  an  excessive  number 
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Table  15.  Incomplete  family  data  collected  from  female  blue  grouse 
and  their  offspring  on  Vancouver  Island,  1969  to  1971. 

The  allele  NgF  is  combined  with  NgS  to  avoid  small  numbers. 


Year 

Mother’s 

Genotype 

Offspring  Genotype 

NgS/NgS 

NgS/NgM 

NgM/NgM 

1969 

NgS/NgS 

3 

3 

NgS/NgM 

4 

4 

7 

NgM/NgM 

- 

8 

25 

1970 

NgS/NgS 

1 

7 

_  * 

NgS/NgM 

3 

19 

27 

NgM/NgM 

-  * 

61 

88 

1971 

NgS/NgS 

3 

7 

_  * 

NgS/NgM 

9 

42 

58 

NgM/NgM 

- 

41 

195 

Totals 

NgS/NgS 

7 

17 

- 

NgS/NgM 

16 

65 

92 

NgM/NgM 

— 

110 

308 

*The  single  offspring  which  belonged  to  each  of  these  cells  was  not 
compatible  with  genetic  theory.  These  individuals  have  been  excluded 
from  all  further  analyses. 
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Table  16.  Distribution  of  the  number  of  young  scored  from  each  mother. 

These  are  the  data  used  to  correctly  calculate  heterogeneity 
X2  values  based  on  a  formula  given  by  Cooper  (1968) .  NgF=NgS. 


Year 

Mother ' s 
Genotype 

Number  of  young 

typed 

1 

2 

3 

4 

5 

6 

^7 

1969 

NgS/NgS 

2 

0 

0 

2 

0 

0 

0 

NgS/NgM 

10 

1 

2 

1 

0 

0 

0 

NgM/NgM 

14 

9 

1 

1 

0 

0 

0 

1970 

NgS/NgS 

1 

1 

0 

0 

1 

0 

0 

NgS/NgM 

8 

2 

2 

1 

4 

0 

1 

NgM/NgM 

17 

11 

6 

5 

8 

1 

3 

1971 

NgS/NgS 

1 

1 

0 

2 

0 

0 

0 

NgS/NgM 

13 

5 

7 

8 

3 

3 

0 

NgM/NgM 

31 

17 

16 

8 

3 

5 

6 

Totals 

NgS/NgS 

4 

2 

0 

3 

1 

0 

0 

NgS/NgM 

31 

8 

11 

10 

7 

3 

1 

NgM/NgM 

62 

37 

23 

13 

11 

6 

9 

■ 


53 


Tabic  17.  Observed  and  expected  genotypic  frequencies  of  mothers  and 

offspring.  Expectations  were  based  on  Levene's  exact  formulae 
for  Hardy- Weinberg  proportions  (Dobzhansky  and  Levene  1948). 

N  is  the  sample  size,  ne  the  effective  number  of  alleles 
(Crow  and  Kimura  1970),  p  the  frequency  of  NgM,  and  probability 
is  the  exact  probability  of  observing  this  distribution  by  chance 
alone.  NgF=NgS. 


Year 

Group 

N 

NgS/NgS 

NgS/NgM 

NgM/NgM 

P 

ne 

probabil 

1969 

Adults 

Obs . 
Exp . 

41 

0.073 

0.055 

0.341 

0.378 

0.585 

0.567 

0.756 

1.58 

0.52 

Chicks 

Obs . 
Exp . 

54 

0.130 

0.068 

0.278 

0.400 

0.593 

0.531 

0.731 

1.65 

0.02 

1970 

Adults 

Obs . 
Exp. 

72 

0.042 

0.026 

0.250 

0.282 

0.708 

0.692 

0.833 

1.38 

0.32 

Chicks 

Obs . 
Exp. 

206 

0.019 

0.052 

0.422 

0.357 

0.558 

0.591 

0.769 

1.55 

0.001 

1971 

Adults 

Obs . 
Exp. 

129 

0.031 

0.032 

0.302 

0.300 

0.667 

0.668 

0.818 

1.42 

0.94 

Chicks 

Totals : 

Obs . 
Exp . 

355 

0.034 

0.025 

0.254 

0.270 

0.713 

0.704 

0.839 

1.37 

0.24 

Adults 

Obs. 

Exp. 

242 

0.041 

0.035 

0.293 

0.307 

0.665 

0.659 

0.812 

1.44 

0.50 

Chicks 

Obs . 
Exp. 

615 

0.037 

0.037 

0.312 

0.313 

0.650 

0.650 

0.807 

1.45 

0.97 

' 

1 
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Table  18.  Results  of  x2  goodness-of-fit  tests  for  the  1:1  ratio  at 
the  Ng  locus  in  blue  grouse.  Ratio  tested: 

°AB-AB:0AB.A+0AB-B  =  1:1* 


Year 

Observed  Numbers 

X2 

1969 

4:11 

2.4 

0.12 

1970 

19:30 

2.0 

0.15 

1971 

42:67 

5.3 

0.02 

Total 

65:108 

10.2 

0.001 
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of  homozygous  offspring.  In  1971,  and  in  the  totals  for  years  combined, 
these  differences  were  significant.  The  excess  of  homozygotes  was  due  to 
zygotic  selection  against  heterozygous  offspring  from  heterozygous  mothers. 

Gametic  gene  ratios  were  tested  and  results  are  presented  in  Table  19. 
These  tests  were  testing  the  gametic  ratios  of  male  parents.  There  were  no 
significant  departures  from  expectations  in  1969  or  1971.  In  1970  rare 
homozygotes  from  heterozygous  parents  were  significantly  different.  In  all 
years  heterozygotes  tended  to  produce  too  few  rare  homozygotes,  resulting 
in  the  totals  being  significantly  different  from  expectations.  That  this 
was  not  consistent  through  all  genotypes  was  shown  by  significant  hetero¬ 
geneity  in  the  gametic  gene  ratios  for  the  totals. 

Based  on  these  results,  departures  from  expectations  in  both  the  1:1 
tests  and  the  gametic  gene  ratio  tests  were  caused  by  selection  against 
heterozygous  and,  perhaps,  rare  homozygous  offspring  of  heterozygous  females. 
Segregation  distortion  can  be  eliminated  as  a  possibility,  since  it  can  be 
shown  that  segregation  distortion  will  not  alter  the  expectations  in  the 
1:1  tests. 

Analogous  tests  to  those  above  can  be  done  on  the  columns  and  diagonals 
of  Table  14.  Using  these  ratios,  the  following  should  hold: 


A-AB  _  °B-AB  _  1 

0  1 
AB-A  AB-B 


0 


0 


AB-AB 


and 


Also,  offspring  from  AA,  AB  and  BB  parents  should  be  in  Hardy-Weinberg 
equilibrium  with  respect  to  each  other  (last  column  of  Table  14).  These 
tests  measure  not  the  viabilities  of  offspring  but  the  reproductive  output 


of  mothers. 
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Table  19.  x2  tests  for  heterogeneity  of  male  gametic  gene  ratios  at 
the  Ng  locus  in  blue  grouse. 


Year  Ratio  Tested 


Observed  Numbers  x2  df  p 


Heterogeneity 
corrected 
for  multiparity 


2  df 


1969  (1) 

°A. A : °A. AB 

3:3 

1.1 

1 

0.3 

(2) 

°AB.A:0AB.B 

4:7 

0.2 

1 

0.7 

(3) 

°B.AB:0B.B 

8:25 

0.5 

1 

0.5 

(4) 

Heterogeneity 

1.8 

2 

0.4 

1.5  p  $  0.47 

1970  (1) 

1:7 

2.1 

1 

0.15 

(2) 

3:27 

9.9 

1 

0.002 

(3) 

61:88 

2.5 

1 

0.12 

(4) 

- 

14.4 

2 

0.001 

9.3  p  $  0.01 

1971  (1) 

3:7 

1.0 

1 

0.3 

(2) 

9:58 

0.6 

1 

0.4 

(3) 

41:195 

0.0 

1 

0.9 

(4) 

- 

1.7 

2 

0.4 

1.3  p  $  0.52 

Total  (1) 

7:17 

0.3 

1 

0.6 

(2) 

16:92 

5.7 

1 

0.02 

(3) 

110:308 

1.1 

1 

0.3 

(4) 

- 

7.1 

2 

0.03 

5.7  p  tj  0.06 
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Results  from  the  auxiliary  1:1  ratio  tests  are  presented  in  Table  20. 
The  first  1:1  ratio  test  showed  a  consistent  deficiency  of  heterozygous 
offspring  produced  by  heterozygous  parents.  All  cases  with  significant 
departures  from  the  expected  ratio  were  a  result  of  too  few  heterozygous 
offspring  being  produced  by  heterozygotes.  This  result  effectively  elimin¬ 
ates  misclassification  as  the  cause  of  the  earlier  heterozygote  deficiency 
and  also  shows  that  selection  was  occurring  only  against  heterozygous  off¬ 
spring  from  heterozygous  mothers.  There  was  only  one  other  significant 
departure  in  the  1:1  ratio  tests.  In  1970,  heterozygous  females  produced 
too  few  common  homozygous  offspring  when  compared  to  the  production  of 
heterozygotes  by  the  common  homozygotes. 

Results  of  the  auxiliary  gametic  gene  ratio  tests  showed  no  consistent 
trends  (Table  21).  Among  the  totals,  rare  homozygotes  produced  fewer 
heterozygous  offspring  than  would  be  expected  when  compared  to  the  output 
of  the  more  frequent  homozygote. 

Summarizing,  the  three  bands  found  on  starch  gels  after  electrophoresis 
have  characteristics  similar  to  other  analogous  systems,  suggesting  the 
bands  were  gene  products  of  three  codominant,  autosomal  alleles.  Incomplete 
family  data  supported  this  hypothesis  and  also  showed  selection  acting 
against  heterozygous  offspring  from  heterozygous  mothers.  These  results 
validate  the  white  bands  as  a  genetic  marker. 

Hardy-Weinberg  Calculations  and  Distribution 

Hardy-Weinberg  Calculations .  The  Hardy-Weinberg  model  relates  observed 
frequency  of  alleles  to  expected  genotype  frequencies  and  makes  comparisons 
of  observed  results  with  expectations  possible.  The  model  is  particularly 
useful  in  the  empirical  study  of  population  genetics  and  is  the  basis  for 
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Table  20. 


Results  of  x2  goodness-of-fit  tests  on  auxiliary  1:1  ratios 
at  the  Ng  locus  in  blue  grouse. 


Year 

Ratio  Tested 

Observed  Numbers 

X2 

P* 

1969 

t1)  °AB. AB  *  °A .AB+0B.AB 

4:11 

2.4 

0.1 

°A.AB • °AB. A 

3:4 

0.0 

1.0 

°B.AB:0AB.B 

8:7 

0.0 

1.0 

1970 

(1) 

19:68 

26.5 

10"4 

(2) 

7:3 

0.9 

0.3 

(3) 

61:27 

12.4 

10"4 

1971 

(1) 

42:48 

0.3 

0.6 

(2) 

7:9 

0.1 

0.8 

(3) 

41:58 

2.6 

0.1 

Total 

CD 

65:127 

19.4 

10"4 

(2) 

17:16 

0.0 

1.0 

(3) 

110:92 

1.4 

0.2 
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Table  21. 


X2  tests  for  heterogeneity  of  auxiliary  gametic  gene  ratios 
at  the  Ng  locus  in  blue  grouse. 


Year 

Ratio  Tested 

Observed  Numbers 

X2 

P* 

1969 

(1) 

°A.A:0AB.A 

3:4 

0.5 

0.5 

(2) 

°A.AB :0B.AB 

3:8 

0.0 

0.8 

(3) 

°AB.B:0B.B 

7:25 

1.1 

0.3 

(4) 

Heterogeneity 

- 

1.2 

0.5 

1970  (1) 

1:3 

0.2 

0.7 

(2) 

7:61 

21.8 

io-4 

(3) 

27:88 

6.9 

0.001 

(4) 

- 

5.4 

0.07 

1971 

(1) 

3:9 

0.5 

0.5 

(2) 

7:41 

0.2 

0.7 

(3) 

58:195 

6.0 

0.01 

(4) 

- 

1.9 

0.4 

Total 

(1) 

7:16 

0.5 

0.5 

(2) 

17:110 

9.1 

0.003 

(3) 

92:308 

0.3 

0.6 

(4) 

7.0 

0.03 
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the  development  of  most  theoretical  population  genetics.  It  is  also  one 
of  the  most  useful  ways  to  reduce  data  and  search  for  patterns  of  deviations 
from  expectations.  Observed  and  expected  allelic  and  phenotypic  frequencies 
for  the  Ng  locus  were  calculated  for  all  sex  and  age  groups  for  all  inten¬ 
sive  study  areas  in  all  years  (Tables  22,  23,  and  24).  There  were  few 
deviations  from  Hardy-Weinberg  expectations  (using  Levene's  exact  expecta- 


4”  i  ■»“'  r—  -f-  ^  r*  *v> 

c  j-iJiiJ  16a  Oil. 


ill  samples  [ Dobzhansky  and  Lcvcne  1948,  Crow  and  Kimura  1970]) 


among  either  adults  or  yearlings.  However,  all  significant  or  near-signif icant 
deviations  were  in  the  direction  of  too  few  heterozygotes. 

Often  the  actual  number  of  alleles  was  greater  than  the  effective  number, 
nfi  (Crow  and  Kimura  1970).  ne  has  also  been  presented  in  these  tables  and 
was  calculated  using  the  relationship 

ne  =  1 


where  p^  is  the  frequency  of  the  i-^  allele  (Crow  and  Kimura  1970).  In 
this  system  ne  had  limits  of  1  and  3.  Most  values  were  between  1.5  and  1.7, 
making  the  effective  number  of  alleles  about  50%  of  the  actual  number. 
Frequency  of  alleles  between  populations  showed  no  consistent  differences. 

The  frequency  of  the  most  common  allele,  NgM>  ranged  from  a  low  of  0.636  to 
a  high  of  1  (in  a  small  sample  of  4)  but  nearly  all  values  clustered  around 
0.8.  NgS  varies  from  0.1  to  0.15  (Figure  8). 

Spatial  and  Temporal  Distribution.  The  same  three  alleles  were  found 
in  approximately  the  same  frequencies  in  populations  of  blue  grouse  on 
Vancouver  Island  from  40  km  south,  100  km  southeast,  30  km  north,  70  km 
northwest,  of  the  Ash  River  Valley,  and  from  mainland  British  Columbia,  110 
km  northeast  of  the  Ash  River  Valley  (Figure  9).  There  was  no  heterogeneity 
in  these  allele  frequencies  between  areas.  The  widespread  occurrence  of  this 
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Table  23.  Genetic  parameters  at  the  Ng  locus  in  blue  grouse  on  area  107,  1968  to  1971.  (Symbols  as  in 
Table  22). 
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Table  24.  Genetic  parameters  at  the  Ng  locus  in  blue  grouse  for  sex  and  age  groups  sampled  on  area  108e, 
1969  to  1971.  (Symbols  as  in  Table  22). 
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Figure  8. 


Frequency  of  the  two  most  common  alleles,  NgM  and  NgSj 
in  blue  grouse  on  the  three  main  study  areas,  1968  to 


1971. 
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Figure  9.  Frequency  histograms  of  the  three  alleles,  NgF  3  NgM  and 
NgSj  in  blue  grouse  from  populations  on  Vancouver  Island 
and  from  the  adjacent  mainland  of  British  Columbia.  The 
three  bars  in  the  histograms  represent,  from  left  to  right 
respectively,  NgS3  NgM  and  NgF.  Sample  sizes  are  given 
under  the  histogram. 

Sample  locations:  1.  Copper  Canyon  (on  Boulder  Creek); 

2.  Cous  Creek;  3.  Stirling  Arm  of  Sproat  Lake;  4.  Taylor 
River  (Taylor  River  Burn);  5.  Ash  River  Valley;  6.  Comox 
Bum  near  Brown's  River;  7.  Middle  Quinsam  Lake;  and 
8.  Pasha  Lake  near  Powell  River. 


100  KM 


HISTOGRAM 

SCALE 


66 


polymorphism  and  the  relationship  of  allele  frequencies  suggest  that  the 
system  was  balanced  (Ford  1965a). 

There  were  some  differences  in  allele  and  genotype  frequencies  between 
areas,  years,  and  sex-age  groupings  in  the  Ash  River  Valley  (Tables  22,  23, 
and  24).  In  1971,  104  was  significantly  different  from  the  other  two  areas 
(X2[2]=8.5,  p<0.01)  because  the  frequency  of  NgM  in  104  was  significantly 
lower  than  in  either  107  (x2[l]=4.17,  p*0.05)  or  108e  (x2[l]=5.9,  p^O.Ol). 

For  areas  combined,  in  1969  there  was  significant  heterogeneity  (x2 [3] =1 1 . 36 , 
p^O.OOl),  because  of  a  high  frequency  of  NgM  among  yearling  males.  Summing 
all  data  for  areas  in  a  given  year  and  comparing  between  years  showed  no 
differences  (x2[3]=6.3,  p^0.093).  All  other  comparisons  showed  no  differ¬ 
ences. 

Similar  calculations  to  those  above  were  done  for  juveniles  (Table  25). 
Since  1968,  juveniles  on  104  had  a  lower  frequency  of  NgM  than  juveniles  on 
either  107  or  108e.  The  frequency  of  NgM  among  juveniles  increased  steadily 
in  108e.  Phenotypic  ratios  among  juveniles  in  1970  (x2[2]=20.3,  p^O.OOOl) 
and  1971  (x2[2]=14.7,  p^0.0006)  were  significantly  different.  The  differ¬ 
ence  was  caused  by  a  low  frequency  of  NgM  among  juveniles  on  104  and  a  high 
frequency  on  108e.  Comparisons  of  phenotypic  ratios  within  areas  but 
between  years  showed  significant  heterogeneity  (x2 [3] =10 . 99 ,  p^0.05)  (104); 
(x2[3]=9.17,  pS0.03)  (107);  (x2[2]=5.9,  pS0.05)  (108e).  In  104  the  hetero¬ 
geneity  was  a  result  of  a  significant  deficiency  of  heterozygotes  in  1971 
(X2 [1]=6.32,  p^O.Ol);  in  area  107  this  was  caused  by  a  low  frequency  of  NgM 
in  1968  (x2[1]=4.69,  p£0.03);  and  in  area  108e  this  resulted  from  a  low 
frequency  of  NgM  in  1969  (x2[l]=4.48,  p£0.03). 

Various  genetic  parameters  for  adults  and  yearlings  combined  were  com¬ 
pared  to  density  of  birds  to  see  if  there  was  a  relationship  between  density 
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and  genetics  on  a  general  level.  There  was  no  correlation  between  adult 
male  density  and  frequency  of  heterozygotes  (r=-0.307),  effective  number  of 
alleles  (r=-0.0442),  0/E  ratio  (observed  number  of  heterozygotes/expected 
number)  (r=-0.282)  or  frequency  of  NgM  (r=-0.2S)  (all  with  11  pairs). 

Summarizing,  nearly  all  Hardy- Weinberg  calculations  showed  no  signi¬ 
ficant  deviations  from  expectations.  The  frequency  of  the  three  alleles 
at  the  Ng  locus  were  widespread  geographically  and  stable  temporally.  There 
were  no  correlations  between  genetic  parameters  and  density.  Thus,  the  Ng 
locus  is  not  only  a  genetic  marker,  it  is  a  stable  marker.  Potentially, 
at  least,  this  increases  the  generality  of  genetic  statements  concerning 
demographic  events. 

Components  of  Fitness 

Various  components  of  fitness  were  analyzed  in  two  Qegmeptc •  ^pro¬ 
duction  and  survival. 

Reproduction.  Fitness  for  breeding  was  analyzed  in  two  ways:  by 
comparing  success  of  different  genotypes  at  completing  nesting  (only  50-80% 
of  females  were  usually  successful  at  nesting),  and  by  comparing  brood  size 
by  maternal  genotype. 

First,  comparisons  of  genotypic  frequencies  of  females  with  broods  to 
genotypic  frequencies  of  those  without  broods  (Table  26)  showed  no  trends 
or  significant  differences  between  successful  and  unsuccessful  breeders; 
(1968--x2[2] =0.996,  psO.617;  1969--X2 [2] =1 . 03 ,  p$0.59;  1970-~X2[2]= 

1.40,  p^0.496;  197.1--X2  [2]  =1 . 09 ,  p$0.58).  Nor  were  there  any  differences 

when  adults  were  considered  separately  from  yearlings  or  when  all  years  were 
combined.  Thus,  success  at  nesting  was  not  selective,  at  least  at  the  tig 


locus . 


, 
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Table  26.  Number  of  female  blue  grouse  of  each  Ng  genotype  with  and 
without  a  brood  in  all  years.  For  computations,  NgF=NgS. 


Year 

Age 

Status 

NgS/NgS 

NgS/NgM 

NgM/NgM 

1968 

Adult 

with 

4 

14 

25 

without 

1 

2 

2 

Yearling 

with 

1 

10 

15 

without 

0 

3 

12 

1969 

Adult 

with 

3 

18 

30 

without 

2 

8 

16 

Yearling 

with 

5 

17 

26 

without 

2 

8 

21 

1970 

Adult 

with 

6 

23 

57 

without 

2 

6 

13 

Yearl ing 

with 

1 

15 

36 

without 

2 

3 

13 

1971 

Adult 

with 

3 

23 

68 

without 

1 

4 

7 

Yearling 

with 

3 

31 

64 

without 

2 

8 

23 

Total 

Adult 

with 

16 

78 

180 

without 

6 

20 

38 

Yearling 

with 

10 

73 

141 

with out 

6 

22 

69 

70 


Second,  comparisons  of  average  brood  size  by  genotype  and  year  were 
made  using  a  two-factor  analysis  of  variance  with  unequal  and  dispropor¬ 
tionate  sample  sizes  by  fitting  constants  (Steel  and  Torrie  1960,  p.257f.). 
Adjusted  average  brood  size  by  age  and  genotype  of  mother  are  presented  in 
Table  27  (adjustments  made  by  fitting  constants).  Interaction  was  insigni¬ 
ficant.  There  were  no  significant  differences  either  between  years  or 
between  genotypes  within  age-classes  (Table  28).  Thus,  the  number  of  chicks 
raised  to  fledgling  stage  was  independent  of  maternal  genotype  at  the  Ng 
locus . 

Survival.  Overwinter  survival  of  juveniles  and  older  birds  was 
analyzed  in  the  following  manner: 

For  juveniles,  the  observed  ratio  of  each  genotype  in  year  X  was  used 
as  the  expectation  among  yearlings  in  year  X+l.  This  assumed  that  yearlings 
in  year  X+l  came  from  a  population  of  juveniles  which  was  statistically  like 
those  sampled  in  year  X.  It  was  not  possible  to  follow  survival  of  marked 
juveniles  since  few  grouse  marked  as  juveniles  were  found  again.  For  older 
birds,  overwinter  survival  of  each  genotype  was  calculated  by  following 
survival  of  tagged  yearlings  and  adults.  This  assumed  that  I  found  marked 
birds  in  the  following  year  at  random  with  respect  to  genotype.  Since  a 
high  percentage  of  birds  was  found  each  year,  this  was  probably  not  invalid. 

Fitness  (W)  values  were  calculated  by  the  formula: 


0.  . 

f 

_il  x 

1 

E.  . 

MAX  0 . . 

il 

E.  . 

13 

where  Ch  .  and  E^_.  are  the  observed  and  expected  proportions,  respectively, 
among  genotype  ij.  The  term  in  brackets  scales  all  fitness  values  to  1  and 
below.  For  juveniles,  the  significance  of  these  values  was  tested  using  x2 
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Table  27.  Adjusted  average  brood  size  of  blue  grouse  by  genotype  and 
year  for  adults  and  yearlings  considered  separately.  The 
adjusted  means  were  calculated  by  the  method  of  fitting 
constants  (Steel  and  Torrie  1960,  p.257f.)#  For  computation 
NgF=NgS .  These  adjustments  were  necessary  for  the  analysis 
of  variance  presented  in  Table  28. 


Age 

Year 

NgS/NgS 

Genotype 

NgS/NgM 

NgM/NgM 

Yearly 

Means 

Adult 

1968 

2.99 

3.32 

3.70 

3.34 

1969 

3.16 

3.48 

3.86 

3.50 

1970 

4.29 

4.62 

5.00 

4.64 

1971 

4.10 

4.40 

4.80 

4.40 

Genotype 

Averages 

3.63 

3.96 

4.33 

3.98 

Yearling 

1968 

3.39 

3.73 

3.44 

3.52 

1969 

2.84 

3.18 

2.89 

2.97 

1970 

3.52 

3.85 

3.56 

3.64 

1971 

3.54 

3.88 

3.58 

3.67 

Genotype 

Averages 

3.32 

3.66 

3.37 

3.45 

, 
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Table  28.  Results  of  two-way  analysis  of  variance  with  unequal,  dis¬ 
proportionate  sample  sizes  on  brood  sizes  and  genotypes  from 
Table  27.  The  method  used  is  that  of  fitting  constants. 
Interaction  is  negligible. 


Age 

Source 

df 

SS 

MS 

F 

Adult 

Genotype  and  Year 

11 

135.70 

12.34 

0.52 

Years  adjusted  for 
genotype 

3 

79.99 

26.66 

1.13 

Genotypes  adjusted 
for  years 

2 

15.14 

7.57 

0.32 

Error 

267 

6513.30 

23.59 

Yearling 

Genotype  and  Year 

11 

31.10 

2.83 

0.18 

Years  adjusted  for 
genotype 

3 

13.95 

4.65 

0.29 

Genotypes  adjusted 
for  years 

2 

3.75 

1.87 

0.12 

188  3018.90  16.06 


Error 


' 
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goodness-of-fit,  while  for  yearlings  and  adults,  significance  was  tested  with 
contingency  x2- 

There  were  no  significant  differences  in  overwinter  survival  among 
juvenile  genotypes  for  1969  to  1970  or  for  1970  to  1971,  but  there  was  a 
significant  difference  in  overwinter  survival  for  juveniles  in  1968  to  1969 
with  NgM/NgM  homozygotes,  the  most  common  genotype,  also  the  most  fit  (Table 
29).  Thus,  juveniles  of  the  most  common  genotype  survived  best  over  one 
winter. 

Among  yearlings  and  adults,  rare  homozygotes,  NgS/NgS3  survived  poorest 
and  heterozygotes  NgS/NgM  survived  best  over  two  winters  (1968  to  1969  and 
1970  to  1971)  (Table  30).  The  winters  1968  to  1969  and  1970  to  1971  were 
severe,  with  the  former  being  one  of  the  hardest  on  record.  In  contrast, 
the  winter  1969  to  1970  was  mild.  In  these  two  severe  winters,  survival 
depended  on  Ng  genotype  but  in  the  mild  winter  survival  was  random  with 
respect  to  the  Ng  locus. 

Summarizing,  selection  did  not  operate  on  reproductive  events,  at  least 
at  the  Ng  locus,  but  selection  occurred  in  overwinter  survival  of  both 
juveniles  and  older  birds.  Juveniles  of  NgM/NgM  genotype  survived  best  over 
one  winter.  In  older  grouse,  heterozygotes,  NgS/NgM3  had  a  selective  advan¬ 
tage  with  respect  to  overwinter  survival,  especially  in  years  of  severe 
winter  weather. 

Genetic  Changes  During  Colonization 

Three  methods  were  used  to  examine  genetic  structure  of  colonizing 
populations :  allelic  frequencies,  frequency  of  heterozygotes  and  an 
observed/expected  (0/E)  ratio  of  heterozygotes.  If  populations  are  in  Hardy- 
Weinberg  equilibrium,  allelic  frequencies  may  give  the  most  pertinent  genetic 
information.  Frequency  of  heterozygotes  was  calculated  by  summing  the 
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Table  29.  Overwinter  survival  of  juvenile  blue  grouse.  Fitness  was 

based  on  observed  and  expected  frequencies  of  each  genotype 
among  yearlings.  Expectations  were  calculated  from  the 
genotype  frequencies  found  among  juveniles  the  previous  summer. 
X2  goodness-of-fit  tests  were  calculated  (2  df)  on  the  observed 
and  expected  numbers.  NgF=NgS  for  computation. 


Year 

N 

NgS/NgS 

NgS/NgM 

NgM/NgM 

x2 

P* 

1969 

Obs . 

116 

0.052 

0.267 

0.680 

13.30 

0.001 

Exp. 

0.073 

0.415 

0.512 

W+ 

0.536 

0.484 

1 

1970 

Obs . 

97 

0.041 

0.258 

0.701 

3.74 

0.15 

Exp. 

0.093 

0.278 

0.630 

Wt 

0.396 

0.834 

1 

1971 

Obs . 

198 

0.04 

0.313 

0.646 

0.95 

0.62 

Exp. 

0.031 

0.338 

0.631 

wt 

1 

0.718 

0.793 

Total 

Obs . 

411 

0.044 

0.287 

0.670 

8.98 

0.01 

Exp. 

0.057 

0.346 

0.597 

Wt 

0.683 

0.740 

1 

t  -  Overwinter  fitness 
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Table  30.  Overwinter  survival  of  yearling  and  adult  blue  grouse. 

Observed  numbers  of  each  genotype  which  survive  and  die 
over  winter;  these  counts  are  based  on  individually  tagged 
birds.  NgF=NgS  for  computation. 


Period 

NgS/NgS 

NgS/NgM 

NgM/NgM 

X2 

P* 

1968-69 

Sur . 

2 

25 

31 

6.99 

0.03 

Die 

6 

14 

43 

Wt 

0.39 

1 

0.65 

1969-70 

Sur. 

7 

47 

87 

0.01 

1.00 

Dae 

6 

40 

76 

Wt 

1 

1 

0.99 

1970-71 

Sur. 

3 

45 

93 

7.02 

0.03 

Die 

13 

40 

91 

Wt 

0.35 

1 

0.95 

Total 

Sur. 

12 

117 

211 

7.00 

0.03 

Die 

25 

94 

210 

Wt 

0.58 

1 

0.90 

t  -  Overwinter  fitness 
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observed  frequencies  of  all  heterozygous  genotypes  within  each  habitat  type. 
Finally,  0/E  ratios  for  each  habitat  type  were  calculated.  The  0/E  ratio 
of  heterozygotes  is  a  measure  of  departure  from  Hardy- Weinberg  expectations. 
If  the  ratio  is  less  than  1, there  is  an  overall  deficiency  of  heterozygotes 
and  if  the  ratio  is  greater  than  1, there  is  a  heterozygous  excess.  The 
null  hypothesis  being  tested,  in  each  case,  was  that  there  were  no  changes 
in  genetic  parameters  associated  with  age  of  habitat  or  population  density. 

Genetic  statistics  on  composition  of  populations  on  different  age 
habitats  since  1968  are  presented  in  Table  31.  There  were  no  significant 
differences  between  sexes,  ages,  or  areas;  all  data  were  combined  for  these 
comparisons.  For  years  combined,  there  was  a  high  frequency  of  NgM  allele 
(0.9)  in  birds  on  one-year-old  habitat,  but  the  frequency  of  this  allele 
decreased  to  about  0.8  by  age-class  2  and  did  not  fluctuate  thereafter. 

Among  totals,  frequency  of  NgM  among  birds  on  age-class  1  was  significantly 
different  from  all  other  age  classes  (Table  32)  but  none  of  the  other  age- 
class  comparisons  were  different  from  each  other. 

As  mentioned  in  the  discussion  on  the  density  of  birds  on  habitats  of 
various  ages,  1971  seemed  to  be  a  year  of  large  numbers  of  birds  moving  onto 
one-  and  two-year-old  habitat.  The  surplus  of  birds  available  to  colonize 
seemed  to  be  larger  in  1971  than  in  other  years.  There  was,  likewise,  in 
1971,  a  significantly  higher  number  of  birds  carrying  the  NgM  allele  on  age- 
class  1  habitat  than  on  other  habitat  types.  Of  25  birds  captured  on  one- 
year-old  habitat  in  1971,  only  one  was  not  an  NgM/NgM  homozygote,  being  an 
NgM/NgS  heterozygote.  Thus,  I  conclude  that  grouse  on  one-year-old  habitat 
were  more  likely  to  be  carriers  of  the  most  common  allele  than  are  grouse 
inhabiting  older  age  habitats. 

Next,  the  observed  frequency  of  heterozygotes  was  plotted  against  age  of 
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Total  N  53  98  131  143  175  115  106  251 

q  0.9  0.8  0.81  0.80  0.81  0.79  0.80  0.80 
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Table  32.  Probability  table  for  comparisons  of  the  number  of  NgM  alleles 
vs.  NgF+NgS  alleles  occurring  in  blue  grouse  on  different  age 
habitats.  When  age-class  1  was  compared  to  other  areas, 
Fisher’s  Exact  Test  was  used  (due  to  small  sample  sizes).  For 
all  other  comparisons,  a  2  x  2  contingency  x2  was  calculated. 


Age  of 

Habitat 

Age 

of  Habitat 

1 

2 

3 

4 

5 

6 

2 

0.02 

3 

0.05 

0.61 

4 

0.03 

0.767 

0.89 

5 

0.06 

0.46 

0.91 

0.70 

6 

0.02 

0.97 

0.70 

0.878 

0.53 

7 

0.04 

0.77 

0.93 

0.93 

0.76 

0.88 

*8 

0.02 

0.74 

0.84 

0.95 

0.62 

0.86 

0.96 


, 


8^ 
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habitat  for  all  years  and  for  the  combined  totals  (Figure  10).  Since  1969, 
frequency  of  heterozygotes  had  a  significant  positive  regression  on  age  of 
habitat.  Clearly,  not  only  was  the  frequency  of  alleles  shifting  during 
colonization,  frequency  of  genotypes  was  also  shifting.  The  slopes  of  these 
regression  lines  were  not  significantly  different  from  each  other. 

To  further  test  the  regression  of  heterozygosity  on  age  of  habitat,  I 
recalculated  the  regression  statistics,  excluding  age-class  1  samples,  since 
age-class  1  had  a  higher  frequency  of  the  most  common  allele  which  lowered 
the  frequency  of  heterozygotes.  When  these  calculations  were  done,  there 
were  still  significant  positive  regressions  of  heterozygosity  on  age  of 
habitat  (total  sample:  b=0.02,  p^O.Ol),  suggesting  that  the  changes  in  geno¬ 
type  frequencies  were  more  than  a  manifestation  of  changes  in  allele  fre¬ 
quencies  . 

Shifts  in  frequency  of  heterozygotes  may  have  resulted  from  differential 
survival  of  genotypes,  selective  recruitment  of  different  genotypes  into  the 
population,  or  both.  Thus,  genetic  statistics  for  each  age  of  habitat  were 
calculated  for  adults  and  yearlings  considered  separately  (Table  33). 
Regression  of  the  observed  frequency  of  heterozygotes  on  age  of  habitat  was 
not  significant  among  adults  (Figure  11).  A  similar  calculation  for  year¬ 
lings  gave  significant  regression  of  the  frequency  of  heterozygotes  on  age 
of  habitat.  However,  the  slopes  of  these  two  regression  lines  were  not 
different  from  each  other. 

0/E  ratios  on  age  of  habitat  were  plotted  separately  for  adults  and 
yearlings  (Figure  12).  Among  adults  the  O/E  ratio  and  age  of  habitat  were 
not  related.  For  yearlings,  the  0/E  ratio  had  a  significant  positive  regres¬ 
sion  on  age  of  habitat.  Again,  however,  the  slopes  of  these  regression  lines 
were  not  different  from  each  other.  Among  yearlings,  there  was  a  significant 
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Figure  10.  Regression  of  observed  frequency  of  heterozygotes  at  the 
Ng  locus  in  blue  grouse  on  age  of  habitat,  1968  to  1971. 
*pS0.05 

**p£0.01 
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Table  33.  Statistics  on  genetic  structure  at  the  Ng  locus  of  adults  and  yearlings  separated  on  different 
age  habitats.  Years  were  combined.  (The  symbols  are  the  same  as  in  Table  31.) 
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Figure  11.  Regression  of  the  observed  frequency  of  heterozygotes  at 
the  Ng  locus  in  blue  grouse  on  age  of  habitat  for  adults 
and  yearlings  considered  separately,  all  years  combined. 
**pj?0.01 
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Figure  12.  Regression  of  the  0/E  ratio  (i.e.,  observed  frequency  of 
heterozygotes/expected  number)  at  the  Ng  locus  in  blue 
grouse  on  age  of  habitat  for  adults  and  yearlings  considered 
separately,  all  years  combined. 
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deficiency  of  heterozygotes  in  habitats  1  and  2  and  a  significant  excess 
in  habitat  7. 

Population  density  increased  with  increasing  age  of  habitat  for  seven 
years  after  logging  (Figure  6).  Thus,  both  density  and  heterozygosity  were 
increasing  with  age  of  habitat.  As  expected,  total  heterozygosity  and 
population  density  were  positively  correlated  (r=0.820,  p^O.Ol).  Frequency 
of  heterozygotes  for  adults  and  yearlings  considered  separately  was  plotted 
against  population  density  (Figure  13).  Among  adults  these  two  variables 
were  not  correlated  (r=0.430,  p^O.05)  but  among  yearlings  these  variables 
were  significantly  correlated  (r=0.838,  p£0.05).  When  analogous  plots  were 
made  for  the  0/E  ratio  and  density  among  adults  and  yearlings  considered 
separately,  yearlings,  but  not  adults,  had  a  significant  positive  relation¬ 
ship  (Figure  14).  These  patterns  were  consistent  whether  total  density, 
adult  density,  or  yearling  density  was  correlated  with  observed  numbers  of 
heterozygotes  or  with  the  0/E  ratio  (Table  34).  Whether  changes  in  popula¬ 
tion  density  were  causing  changes  in  allele  and  genotype  frequencies  is  not 
known;  nor  is  the  converse  known.  Perhaps  both  were  caused  by  some  other 
common  factor  such  as  changing  quality  of  vegetation  or  habitat  selection. 
There  can  be  little  question  that  there  were  genetic  shifts  associated  with 
changes  in  abundance  during  colonization. 

The  regression  of  frequency  of  heterozygotes  and  the  0/E  ratio  on  age 
of  habitat  among  yearlings  shows  that  yearlings  were  selectively  recruited 
into  populations  on  different  habitat  types.  The  mechanism  of  selection 
is  not  known,  but  it  seems  unlikely  that  the  birds  were  selecting  some 
environmental  variable  regardless  of  population  density.  Population  density 
consistently  changed  through  the  first  seven  years  after  logging.  These 
results  arc  evidence  in  favor  of  the  view  that  cryptic  genetic  polymorphisms 
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Figure  13.  Correlation,  for  adults  and  yearlings  considered  separately, 
between  the  observed  frequency  of  heterozygotes  and  total 
density  of  blue  grouse.  Circles  are  points  for  adults  and 
triangles  are  for  yearlings. 
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Figure  14.  Correlation,  for  adults  and  yearlings  considered  separately, 
between  the  0/E  ratio  and  total  density  of  blue  grouse. 
Circles  are  points  for  adults  and  triangles  are  for  year¬ 
lings  . 
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Table  34.  Correlation  coefficients  for  comparisons  between  frequency  of 
heterozygotes  and  the  0/E  ratio  compared  to  total t  adult  and 
yearling  density,  all  with  11  pairs  compared. 


Frequency  of  Heterozygotes 

0/E  Ratio 

Adult 

Yearling 

Adult 

Yearling 

Total  Density 

0.426 

0.838** 

0.124 

0.844** 

Adult  Density 

0.533 

0.861** 

0.148 

0.891** 

Yearling  Density 

0.214 

0.693* 

0.248 

0.694* 

*p*0.05 

**p$0.01 
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are  not  selectively  neutral. 

Summarizing,  the  frequency  of  NgM  was  significantly  higher  in  birds 
on  one-year-old  habitat  than  on  other  age  habitats.  The  observed  frequency 
of  heterozygotes  had  a  positive  regression  on  age  of  habitat  among  year¬ 
lings  but  not  among  adults.  The  same  was  true  of  the  0/E  ratio  of  hetero¬ 
zygotes.  Population  density  was  correlated  with  both  heterozygosity  and 
the  0/E  ratio  in  yearlings  but  not  in  adults.  These  results  are  in  agreement 
with  the  original  hypothesis  that  there  would  be  genetic  changes  associated 
with  changes  in  population  density. 
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DISCUSSION 

Colonization  and  Population  Regulation 

This  study  was  designed  to  describe  colonization  by  blue  grouse  of 
newly  logged  regions  (a  period  when  populations  of  grouse  are  rapidly  ex¬ 
panding)  and  to  test  the  hypothesis  that  genetic  structure  of  populations 
would  shift  during  this  episode.  The  hypothesis  stems  directly  from  Chitty's 
(1967)  suggestion  that  genetic  changes  occur  during  changes  in  population 
density. 

A  fundamental  question  in  population  biology  is  why  do  most  populations 
tend  to  remain  stable,  i.e.,  why  are  rates  of  increase  over  long  periods 
near  to  zero?  Generally,  animals  produce  more  offspring  than  necessary  to 
replace  adult  mortality  (Lack  1954,  1966,  and  others).  This  production  is 
a  constant  driving  force  away  from  equilibrium.  Population  regulation  is 
the  process  which  counteracts  this  force  and  promotes  stability  (Chitty 
1960,  Krebs  1964,  Schwerdtefeger  1968). 

Several  theories  have  been  advanced  to  account  for  the  regulation  of 
numbers  of  animals  in  natural  populations.  Lack  (1966)  and  others  suggest 
that  population  regulation  is  brought  about  by  external  factors,  especially 
the  amount  of  food  available  outside  the  breeding  season.  In  their  view, 
animals  are  reproducing  as  rapidly  as  possible  and  population  regulation 
occurs  because  of  factors  beyond  control  of  the  population.  This  has  been 
called  extrinsic  regulation  of  numbers. 

Wynne- Edwards  (1962)  suggests  that  animals  limit  their  reproductive 
output  by  means  of  social  behaviour  and  assessment  of  resources.  This  is  a 
self-imposed  birth  control  system,  geared  to  external  factors  and  involves 
the  concept  of  group  selection. 
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Chitty  (1967),  arguing  in  a  somewhat  similar  vein  but  not  involving 
the  group  selection  of  Wynne- Edwards ,  suggests  that  numbers  are  regulated 
by  genetically  controlled  behaviour,  i.c.,  populations  are  self-regulating 
because  certain  genotypes  have  an  advantage  when  populations  are  sparse  but 
not  when  they  are  dense.  Both  Chitty' s  and  Wynne- Edwards '  hypotheses  have 
been  called  intrinsic  regulation  of  numbers.  However,  both  Chitty  and 
Wynne- Edwards  recognize  the  need  for  external  factors  to  act  as  environmental 


cues 


The  results  of  extrinsic  vs.  intrinsic  control  are  quite  obvious.  If 
a  population  is  regulated  by  extrinsic  mechanisms,  breeding  surpluses  never 
exist.  Thus,  the  colonization  of  new  habitat  should  be  by  drifters  away 
from  established  populations  and  the  populations  which  colonizers  leave 
should  decline.  However,  intrinsic  regulation  allows  for  a  breeding  surplus 
to  be  generated,  since  the  behavioural  interactions  usually  thought  to 
regulate  numbers  take  place  with  the  onset  of  breeding.  Kluiver  and  Tinbergen 
(1960)  demonstrated  a  behavioural  surplus  in  great  tits  (Parus  major)  and, 
using  Chitty' s  hypothesis,  these  surpluses  should  be  genetically  different 
from  non-surplus . 

In  several  more  recent  field  studies  on  the  ecology  and  genetics  of 
small  mammals  (rodents),  attempts  were  made  to  relate  genetic  composition 
to  changes  in  density.  Semeonoff  and  Robertson  (1968)  reported  a  change  in 
the  frequency  of  an  esterase  allele  (Es  locus)  related  with  population  changes 
in  Micro tus  agrestis.  Others  also  reported  genetic  changes  associated  with 
demographic  events.  Canham  (1969),  Gaines,  et  al.  (1971),  Gaines  and  Krebs 
(1971),  Myers  and  Krebs  (1971)  and  Tamarin  and  Krebs  (1971)  all  showed  a 
relationship  between  population  increases  or  decreases  and  genetic  structure 
in  several  species  of  mice. 
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All  studies  outlined  above  were,  concerned  with  changes  in  genetic 
structure  of  populations  with  changes  in  density.  Some  have  shown  that  fit¬ 
ness  and  density  cannot  be  considered  separately,  i.e.,  fitnesses  of  geno¬ 
types  change  with  changes  in  density,  supporting  Kojima’s  (1971)  conclusion 
that  no  single  set  of  fitness  values  can  be  attached  to  a  particular  locus. 

Thus,  Chi tty's  hypothesis,  originally  formulated  to  explain  numbers 
of  small  rodents,  but  broadened  to  include  all  animals  (Chitty  1967),  has 
been  tested  on  only  a  few  populations  of  small  rodents.  Clearly,  the  need 
is  to  go  beyond  this  limited  scope  and  test  the  hypothesis  on  a  wider  variety 
of  organisms. 

In  blue  grouse  the  Ug  locus  is  polymorphic  and  apparently  of  simple 
inheritance.  The  three  alleles  at  this  locus  were  widespread  in  similar 
frequencies  among  birds  on  Vancouver  Island  and  among  birds  on  the  adjacent 
mainland.  Three  lines  of  evidence  suggest  that  this  locus  (or  a  closely 
linked  one)  was  of  some  biological  significance  in  the  life  of  grouse. 

First,  selection  acted  against  heterozygous  offspring  from  heterozygous 
mothers.  While  no  mechanism  is  known  to  account  for  this  phenomenon,  a 
similar  finding  was  reported  for  domestic  cattle  (Bos  taurus)  at  the  Tf 
locus  (Ashton  1965,  Cooper  and  Rendell  1968).  Second,  after  one  year  of  age 
heterozygotes  survived  better  than  homozygotes  during  harsh  winters..  Finally, 
birds  of  different  genotypes  selectively  recruited  onto  habitats  of  various 
ages.  Selective  recruitment  was  correlated  with  population  density,  but 
it  was  not  known  if  the  relationship  was  cause  and  effect. 

The  locus  used  in  this  study  was  chosen  in  a  nearly  random  manner. 

The  same  locus  was  found  in  all  species  of  birds  I  examined  (robins,  Turdus 
migratorius ;  ruffed  grouse,  Bortasa  umbellus;  white- tailed  ptarmigan,  Lagopus 
leucurus;  spruce  grouse,  Canachites  canadensis;  and  sharp- tailed  grouse, 
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Vedi-oecetes  phasianattus )  as  well  as  a  variety  of  mammals  (see  Birds  a  11  et  at. 
1970).  Even  though  nothing  is  known  about  the  biochemical  significance  of 

the  Ng  locus,  -this  does  not  invalidate  it  as  a  genetic  marker.  My  original 
hypothesis  was  stated  in  general  terms  to  cover  changes  in  genetic  structure 
during  increases  in  density.  In  many  respects  the  random  choosing  of  this 
locus  generalizes  the  results. 

In  order  to  test  a  hypothesis  that  genetic  structure  changes  with 
changes  in  density,  two  experimental  approaches  could  be  used  to  induce  rapid 
population  increases  in  normally  stable  populations.  Classically,  the  method 
used  to  induce  population  increases  has  been  the  removal  of  resident  breeders 
(Bendell  and  Elliott  1967,  Harris  1970,  Hensley  and  Cope  1951,  Krebs  1970, 
1971,  Myers  and  Krebs  1971,  Stewart  and  Aldrich  1951,  and  Watson  and  Jenkins 
1968).  Essentially,  removal  experiments  involve  the  creation  of  depopulated 
areas  in  suitable  habitat  and  their  subsequent  repopulation.  A  second  means 
to  induce  rapid  population  increase  is  to  create  suitable  habitat  out  of 
formerly  unsuitable  habitat. 

Colonization  of  unoccupied  habitat  is  a  problem  in  evolutionary  theory 
(Mayr  1963,  1965,  Simberloff  and  Wilson  1969)  and  population  dynamics.  The 
specific  details  concerning  the  mechanisms  of  colonization  have  not  been 
examined  in  many  instances,  but  the  results  of  such  examinations  are.  impor¬ 
tant  to  population  and  evolutionary  biologists  alike.  As  population  biolo¬ 
gists,  the  types  of  questions  we  must  answer  are:  What  is  the  rate  of  move¬ 
ment  onto  unoccupied  regions?  What  is  the  sex  and  age  composition  of 
colonizing  populations?  How  long  is  a  colonizing  episode  for  a  particular 
species?  Are  colonizers  behaviourally  distinct  from  non-colonizers?  Are 
there  genetic  differences  between  colonizers  and  non-colonizers?  Are 
colonizers  surplus  to  other  populations? 
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Density  of  blue  grouse  increased  steadily  for  at  least  seven  years 
following  logging.  Densities  of  adult  males  and  adult  females  increased 


more  rapidly  than  those  of  yearlings.  Yearlings  seemed  to  have  two  stages 
of  colonization.  The  first  was  a  four  year  period  of  stability  at  about 
12  females/km2  and  the  second  was  for  at  least  3  years  at  about  22  females/ 
km  .  Because  populations  of  grouse  in  the  general  vicinity  were  .increasing. 


at  least  ^ ince  19u9 , 


I  cannot  say  whether  colonization  would  always  take 


seven  years.  In  fact,  if  colonization  is  that  period  when  age  structure  of 
populations  is  unstable,  then  colonization  was  over  in  three  to  four  years, 
since  this  was  the  length  of  time  it  took  the  age  structure  to  stabilize. 

It  may  be  difficult  to  distinguish  between  a  colonizing  population  and  one 
that  is  increasing.  In  fact,  the  two  phenomena  are  similar  and  the  mechan¬ 
isms  causing  both  may  be  similar.  Both  involve  the  recruitment  of  young, 
both  have  increasing  densities  and  both  may  involve  genetic  changes. 

The  origin  of  colonizing  populations  was  not  known ,  since  few  grouse 
banded  as  juveniles  were  seen  later.  On  the  basis  of  banding,  it  seems  safe 
to  conclude,  however,  that  most  colonizers  were  produced  outside  the  study 
areas . 

Throughout  the  period  of  population  increase  there  were  systematic 
shifts  in  the  genetic  structure  of  populations.  These  shifts  involved 
changes  in  both  allele  and  genotype  frequencies.  Birds  on  one-year-old 
habitat  had  a  higher  frequency  of  NgM  than  birds  on  all  older  age  habitats 
but  there  were  no  other  changes  in  allele  frequencies.  Heterozygosity  also 
increased  with  increases  in  density.  In  addition  to  these  changes  in  allele 
and  genotype  frequencies , genetic  structure  changed  with  increases  in  density 
from  populations  with  a  heterozygote  deficiency,  to  ones  with  a  balanced 
genetic  structure,  to  ones  with  a  heterozygote  excess.  Thus,  genetic  structur 
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in  these  increasing  populations  has  systematically  changed  from  ones  with 
unbalanced  heterozygous  deficiencies  to  ones  with  unbalanced  heterozygous 
excesses . 

Are  colonizers  surplus  to  other  populations?  The  possibility  exists 
that  colonizers  were  surplus  from  habitat  which  was  no  longer  suitable. 
(Bendell  and  Elliott  [1966]  showed  that  yearlings  do  not  settle  in  dense 
habitat.)  Two  possibilities  exist:  either  most  colonizers  were  surplus 
to  stable  breeding  populations  or  they  were  refugees  from  unsuitable  habitat. 
A  simple  model  can  be  used  to  help  resolve  this  problem.  Assume  that  repro¬ 
duction  was  the  same  per  individual  regardless  of  habitat  type.  (There  is 
no  evidence  from  other  studies  of  blue  grouse  which  would  contradict  this 
assumption.)  If  X  was  the  amount  of  habitat  created  by  logging  each  year, 
and  if  Y  was  the  amount  lost  to  regrowth,  then,  if  X=Y,  production  could  have 
balanced  recruitment  without  an  additional  surplus  generated  by  interaction. 
In  the  decade  1960  to  1969  the  average  area  logged  in  the  Ash  River  Valley 
was  50%  greater  than  in  the  decade  1950  to  1959.  Since  habitat  most  likely 
remains  suitable  for  10  to  15  years  in  this  region,  the  amount  of  new  habitat 
created  in  the  last  few  years  is  greater  than  the  amount  lost  to  regrowth. 
Using  the  solution  given  above,  I  conclude  blue  grouse  colonizing  newly 
logged  regions  were  surplus  to  breeding  populations. 

At  this  point,  this  study  lends  credibility  to  Chitty's  hypothesis, 
but  opens  even  more  far  reaching  problems:  If  a  locus,  chosen  at  random, 
is  associated  with  population  dynamics,  then,  since  at  least  several  thousand 
loci  are  polymorphic  (Lewontin  1967),  several  hundred  or  thousand  polymorphic 
loci  must  also  show  changes  with  changes  in  density.  A  suitable  null  hypo¬ 
thesis  for  any  future  study  on  any  species,  it  seems,  would  be  that  no  poly¬ 
morphic  locus  is  independent  of  demographic  events. 
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One  possible  way  to  resolve  this  problem  might  be  the  following  hypo¬ 
thetical  model:  Individuals  with  the  greatest  degree  of  heterozygosity  have 
the  greatest  competitive  advantage.  Ignoring  the  Ng  locus,  for  the  moment, 
individuals  will  have  a  certain  average  heterozygous  background.  For 
example,  say  an  average  individual  has  200  heterozygous  loci,  excluding  the 
Ng  locus.  Then  individuals  heterozygous  at.  the  Ng  locus  are,  on  average, 
slightly  more  heterozygous  than  homozvgotes  at  this  locus  (201  vs.  200).  If 
overall  heterozygosity  gives  a  competitive  advantage,  then  Ng  heterozygotes, 
on  average,  would  be  competitively  superior  to  Ng  homozygotes. 

Now,  imagine  that  the  quality  of  the  habitat  of  blue  grouse  increases 
for  several  years  following  logging.  Individuals  compete  in  the  highest 
quality  habitat  first,  but  because  of  the  superiority  of  heterozygotes,  the 
highest  quality  habitat  would  be  filled  by  an  excess  of  heterozygotes.  The 
next  best  habitat  would  be  filled  in  a  similar m anner  but  since  there  would 
be  fewer  heterozygotes  available  there  would  be  relatively  fewer  selecting 
this  habitat.  This  process  would  continue  until  all  available  habitat,  was 
filled,  or  until  no  more  birds  were  left.  This  stepwise  settling  of  habitat 
potentially  could  create  a  relationship  between  overall  heterozygosity  and 
age  of  habitat  like  observed.  The  Ng  locus,  then,  becomes  not  necessarily 
important  in  itself  but  a  marker  in  relation  to  a  whole  genetic  background. 

Others  have  shown  that  colonizers  are  more  homozygous  than  non-colonizers. 
As  pointed  out  by  Anderson  (1970) ,  Semeonoeff  and  Robertson  (1968)  found  that 
Microtus  agvestis  inhabiting  temporarily  suitable  habitat  (due  to  periodic 
flooding)  had  a  higher  frequency  of  an  Es  allele  (which  also  meant  there  were 
more  homozygotes)  than  animals  in  more  stable  areas.  Also,  Myers  and  Krebs 
(1971)  found,  most  consistently,  in  M.  ochrogaster  and  M.  pennsylvanicus } 
homozygous  genotypes  dispersed  more  frequently  than  heterozygous  genotypes, 
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especially  during  periods  of  population  increase. 

The  observed  relationship  between  survival  value  and  Ng  genotype  among 
adult  and  yearling  blue  grouse  depended  on  the  severity  of  winter  weather. 
During  two  harsh  winters  heterozygotes  survived  best,  while  during  one  mild 
winter  all  genotypes  survived  equally  well.  Selection  dependent  on  weather 
may  be  a  means  of  causing  population  increases  or  declines.  For  example, 
several  harsh  winters  in  succession  might  change  the  genetic  structure  of 
populations  by  continual  favoring  of  heterozygotes.  If  heterozygotes  are 
competitively  superior  and  able  to  tolerate  crowding,  a  population  increase 
may  follow.  On  the  other  hand,  several  mild  winters  may  mean  that  less 
tolerant  homozygotes  would  become  more  abundant  and  populations  might  decline 

Evolutionary  Implications  of  Experimental  Population  Genetic  Studies 

Harris  (1966),  Lewontin  and  Hubby  (1966),  O'Brien  and  MacIntyre  (1969), 
Prakash  et  al.  (1969)  and  Selander  et  al.  (1969)  have  shown  that  a  large 
proportion  of  the  genome  of  animals  is  polymorphic  with  the  same  sets  of 
alleles  widespread  at  relatively  constant  frequencies.  One  of  the  more 
interesting  outgrowths  of  the  demonstrations  of  large  amounts  of  genetic 
variation  is  the  continuing  discussion  on  the  m aintenance  of  such  variability 
Manwell  and  Baker  (1970),  Prakish  et  al .  and  others  have  argued  that  such 
variability  has  biological  patterns  suggesting  maintenance  by  selection. 
Kimura  (1968,  1969),  Kimura  and  Maruyama  (1971),  and  Kimura  and  Ohta  (1971) 
argue  that  most  polymorphisms  are  neutral  and  that  mutations,  drift  and 
movement  are  enough  to  maintain  the  same  allelic  states  over  wide  areas 
(Kimura  and  Ohta  1971,  Maruyama  1970a, b).  Population  genetics  today  is 
concerned  largely  with  explaining  the  maintenance  of  genetic  variability 
(Wallace  1970).  The  problem  arises  because  if  most  polymorphic  loci  are 
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maintained  by  selection,  and  if  there  are  several  thousand  polymorphic  loci 
(Lewontin  1967),  a  genetic  load  is  placed  on  a  population  which  is  too  great 
to  bear  (Kimura  and  Crow  1964) . 

The  issues  in  this  debate  seem  clear  cut.  Either  most  polymorphisms 
are  maintained  by  selection  or  most  are  not.  If  most  are  not,  then  the 
probability  of  choosing  a  polymorphism  at  random  and  showing  any  differential 
survival  and/or  reproduction  is  vanishingly  small.  Intensive  genetic  and 
demographic  studies,  such  as  this  one,  provide  useful  data  to  this  debate. 
Canham  (1969),  Gaines  and  Krebs  (1971),  Gaines  et  at.  (1971),  Myers  and 
Krebs  (1971),  Semeoneoff  and  Robertson  (1968),  Tamarin  and  Krebs  (1969)  and 
this  study  showed  differential  survival  and/or  reproduction  associated  with 
a  variety  of  relatively  randomly  chosen  cryptic  enzymatic  and  protein  poly¬ 
morphisms.  This  seems,  on  the  surface,  to  be  conclusive  evidence  against 
the  neutral  allele  hypothesis. 

To  circumvent  the  problem  associated  with  genetic  load  arising  from 
non-neutral  loci.  King  (1967),  Milkman  (1967)  and  Sved  et  at.  (1967)  argued 
for  truncation  selection  of  individuals,  not  loci.  Wallace  (1968a, b,  1970) 
proposed  the  concept  of  hard  and  soft  selection  to  account  for  the  observed 
facts . 

Another  way  to  maintain  large  amounts  of  genetic  variability  without 
genetic  load  is  through  multiple-niche  polymorphism.  The  theoretical  develop¬ 
ment  of  this  mechanism  has  been  developed  by  Levene  (1953,  1967)  and  Prout 
(1968).  For  example,  if  animals  of  different  genotypes  select  habitat  at 
random,  but  survive  differentially  according  to  genotype  on  each  habitat,  it 
is  possible  to  have  a  stable  polymorphism  without  heterozygous  advantage. 

In  a  species  such  as  blue  grouse,  adapted  to  living  in  temporarily  suitable 
habitat,  this  might  be  a  means  of  adding  stability  to  a  polymorphism  without 
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increasing  the  genetic  load. 

It  is  likely  that  grouse  are  selecting  habitat  according  to  genotype 
rather  than  being  selected  against  once  they  arrive  onto  a  habitat.  The 
mechanism  of  habitat  selection  by  blue  grouse  is  not  known.  What  is  it 
about  the  Ng  locus  that  makes  it  important  in  this  type  selection?  Is  it 
quality  or  quantity  of  food?  Or  are  there  behavioural  interactions  involved? 
Perhaps  older  habitat  is  most  suitable  and  heterozygotes  are  able  to  preempt 
a  spot  to  live.  Homozygotes,  on  the  other  hand,  may  not  be  able  to  tolerate 
high  density  populations  and  move  to  lower  density  regions. 

This  discussion  illustrates  one  of  the  ways  that  detailed  demographic 
and  genetic  studies  can  aid  in  our  understanding  of  events  occurring  in 
evolution  and  help  resolve  problems  arising  in  theoretical  development  of 
population  biology. 

Validity  of  the  Data 

In  any  field  study,  a  wide  variety  of  uncontrolled  variables  can  be 
disturbingly  annoying  and  this  study  is  not  unique  in  that  respect.  The 
experimental  design  left  something  to  be  desired.  However,  I  had  no  control 
over  logging  patterns  and,  therefore,  no  control  over  the  size,  shape  or 
placement  of  new  habitat.  The  size  of  openings  ranged  from  15  to  over  40 
hectares.  Some  of  the  openings  were  created  in  autumn  or  winter  while  others 
were  created  in  spring  and  summer.  Sometimes  these  new  regions  were  burned 
after  logging;  sometimes  they  were  not.  However,  burning  appeared  to  have 
little  effect  on  density  of  blue  grouse  (Redfield  et  at.  1970).  These  un¬ 
controlled  experimental  variables  probably  made  the  results  less  clear  cut. 
Thus,  the  patterns  that  emerged  were  probably  indicative  of  the  overriding 
pattern  of  colonization. 
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For  purposes  of  analysis  of  colonization,  I  counted  birds  only  once 
in  a  given  year  and  each  bird  was  assigned  to  the  habitat  type  on  which  it 
was  first  found.  A  bird  could  have  been  counted  on  a  given  type  of  habitat 
when  it  belonged  to  another.  For  adult  males,  who  occupy  a  small  territory 
in  spring  and  summer,  this  was  probably  insignificant.  For  other  segments 
of  the  population,  however,  misclassification  was  potential.  I  assumed  that 
any  misclassification  would  be  balanced  by  an  opposite  misclassification. 

Again,  this  might  add  variability  to  the  results  and  obliterate  some  not  so 
obvious  trends.  Thus,  any  relationships  that  emerge  must  be  fairly  indica¬ 
tive  of  what  is  actually  happening. 

Another  source  of  variation  was  the  possibility  of  misclassification 
of  genotypes.  While  there  may  have  been  some  misclassification  of  genotypes 
(both  because  of  clerical  errors,  which  are  impossible  to  eliminate  completely, 
and  also  because  of  leakage,  or  misreading),  of  over  600  chicks  whose  mothers 
were  known,  only  three  were  not  compatible  with  genetic  hypothesis  and  it  is 
possible  to  account  for  these  by  mixing  of  broods.  Also,  most  samples  of 
blood  were  classified  more  than  once.  Even  if  some  misclassification  occurred, 
it  should  have  been  at  random  over  all  the  samples  regardless  of  the  habitat 
type  from  which  the  bird  was  taken. 

One  might  question  the  conclusion  that,  based  on  a  single  locus,  genetic 
structure  changes  during  a  colonizing  episode.  If  this  is  the  only  locus 
which  shows  the  effect  claimed,  then  indeed,  genetic  structure  has  changed. 

If,  however,  other  loci  show  changes  (of  any  nature) ,  then  the  results  are 
even  more  striking.  These  results  force  one  to  say  that  genetic  structure 
has  indeed  changed  during  colonization. 

The  Future 

The  results  of  this  study  are  consistent  with  Chi tty's  genetic  hypothesis, 
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but  they  are  also  consistent  with  other  hypotheses.  The  possibility  that 
birds  selected  features  of  the  habitat  related  to  logging  history  in  relation 
to  genotype  cannot  be  ruled  out.  However,  acceptance  of  this  hypothesis  is 
more  difficult  when  one  considers  that  grouse  select  habitat  in  spring  before 
new  growth  has  started,  and  often  when  snow  is  still  on  the  ground.  If  grouse 
are  selecting  habitat,  they  must  anticipate  the  changes  which  will  occur  over 
the  next  few  months.  A  more  direct  measure  on  which  to  base  selection,  it 
seems,  is  one  at  least  partially  based  on  density  and  kind  of  conspecifics . 

Two  suitable  experimental  approaches  could  test  the  alternatives  outlined 
above.  First,  one  could  continually  remove  resident  populations  for  several 
summers  in  different  habitat  types.  Recruitment  of  new  populations  would  then 
always  be  into  low  density  situations  but  habitat  would  be  changing  between 
years.  Also,  if  care  were  taken  at  the  start,  habitat  types  of  removal  plots 
would  be  different  within  years. 

A  second  experimental  approach  involves  continual  manipulation  of  the 
habitat  in  such  a  manner  as  to  keep  it  in  an  early  stage  of  succession.  For 
example,  one  could  manipulate  by  continually  burning  a  section  of  habitat 
every  year  or  by  continually  cutting  out  all  coniferous  and  deciduous  trees. 

If  genetic  structure  of  recruits  was  determined  by  habitat  and  not  population 
density,  then  one  would  expect  recruits  onto  these  m anipulated  plots  to  be 
genetically  like  early  colonizers.  On  the  other  hand,  if  genetic  structure 
of  recruits  was  determined  in  part  by  population  density,  recruits  onto  these 
manipulated  habitats  would  be  genetically  like  late  colonizers  (provided,  of 
course,  the  populations  did  not  decline  as  a  result  of  habitat  manipulation). 

The  first  experiment  is  rather  practical  in  that  manipulation  of  the 
habitat  is  not  necessary.  Habitat  manipulation  in  regions  where  trees  are 
commercially  valuable  is  difficult  to  justify.  On  the  other  hand,  if  the 
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relationship  between  habitat  and  populations  is  to  be  fully  understood,  we 
must  directly  manipulate  the  habitat  in  a  predetermined  manner.  Thus,  the 
second  experiment  is  every  bit  as  vital  to  these  studies,  but  may  involve 
cooperation  on  a  larger  scale  than  ever  before  in  grouse  research.  This 
cooperation  may  need  to  be  between  biologists,  logging  companies  and  provincial 
governments . 

In  addition  to  these  experimental  approaches,  new  genetic  markers  must 
be  studied.  As  many  more  loci  as  feasible  should  be  studied  in  the  future. 

We  need,  in  short,  a  large  survey  of  the  genome,  as  done  by  Lewontin  and 
Hubby  (1966),  Prakish  et  al.  (1969)  and  Selander  (1970a, b) ,  coupled  with  an 
intensive  population  study.  Only  then  will  the  full  interaction  of  the  genome 
and  density  begin  to  be  revealed. 
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CONCLUSIONS 

What  has  this  research  accomplished  and  revealed? 

1)  Blue  grouse  colonized  newly  logged  regions  rapidly  and  increased 

in  density  for  at  least  seven  years  after  logging.  The  increase  was  brought 
about  by  increases  in  adult  and  yearling  density.  Age  structure  of  colonizing 
populations  shifted  significantly  from  70%  to  45%  yearlings  in  the  first 
three  years  after  logging,  but  stabilized  after  that. 

2)  Reproductive  output  was  high  and  early  colonizing  populations  did 
not  appear  to  be  any  more  successful  at  reproduction  than  late  colonizers. 

3)  A  series  of  white  bands  from  serum  of  grouse  were  identified  with 
starch-gel  electrophoresis.  These  bands  were  inherited  as  codominant  alleles 
at  a  single  autosomal  locus.  Incomplete  family  data  analysis  demonstrated 
that  heterozygous  mothers  produced  too  few  heterozygous  offspring.  No 
mechanism  was  known  which  could  account  for  this,  but  segregation  distortion 
was  ruled  out. 

4)  Birds  of  the  various  genotypes  at  the  Ng  locus  did  not  show  any  dif¬ 
ferential  reproduction  but  survival  was  selective.  First  year  survival 
favored  the  common  homozygote,  NgM/NgM ,  over  one  winter  and  survival  of  older 
birds  favored  heterozygotes,  NgS/NgM,  during  two  severe  winters. 

5)  Frequency  of  heterozygotes  and  the  0/E  ratio  (observed/expected 
heterozygotes)  among  yearlings  had  a  significant  positive  regression  on  age 
of  habitat . 

6)  Correlation  between  population  density  and  both  frequency  of 
heterozygotes  and  the  0/E  ratio  among  yearlings  was  positive  and  significant. 
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Whether  these  shifts  were  caused  by  changes  in  population  density  or  changes 
in  habitat  is  not  known,  but  it  is  difficult  to  imagine  that  the  shifts  were 
occurring  independently  of  population  density.  Suitable  experiments  were 
briefly  outlined  which  could  be  used  to  approach  this  problem. 

7)  The  results  of  this  study,  along  with  those  of  other  field  studies, 
are  not  in  agreement  with  Kimura’s  (1968)  neutral  allele  hypothesis,  since, 
if  Kimura  is  right,  the  probability  of  choosing  a  locus  whose  genotypes  have 
any  differential  significance  to  the  life  of  an  individual  is  small. 

8)  These  results  lend  some  credibility  to  the  hypothesis  concerning 
genetic  changes  associated  with  changes  in  density,  but  at  the  same  time 
raise  some  even  more  perplexing  problems,  not  the  least  of  which  is  the  sug¬ 
gestion  that  several  hundred  or  several  thousand  polymorphic  loci  must  be 
related  with  changes  in  population  density. 
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